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FOREWORD 


' The ma^or computational computer programs described in this manual 
were developed by NASA/Ames Research Center personnel. The Finite-Volume 
Unsteady Blunt-Body program and the Finite-Difference Shock-Capture pro- 
gram were each coded in the CFD language by Ames personnel for use on the 
ILLIAC IV computer. The Finite-Volume Steady Blunt-Body program was 
coded in FORTRAN by Ames personnel and then translated into CFD by the 
authors of this manual. The Coordinate Transformation Interface program 
was also developed by the authors to allow the main space-shuttle programs 
to be run in sequence. All CFD coded programs have been translated back 
into FORTRAN by the CPDX translator for the CDC 7600 at Ames. 4-’ 

The assistance and cooperation of the Technical Monitor, Dr. Harry E. 
Bailey, in providing access to the programs and data needed to accomplish 
the task of combining the required programs to develop a complete space- 
shuttle computation method is gratefully acknowledged. The assistance of 
Dr. Walter Reinhardt, Mr. William Davy and Mr. John Rakich in developing 
the interfaces between the various programs is also gratefully acknowl- 
edged. 
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USERS MANUAL FOR SPACE-SHUTTLE 


COMPUTER PROGRAMS 


By Gary D. Kuhn, Frederxck K. Goodwin 
and Stanley C, Perkins, Jr. 
Nielsen Engineering & Research, Inc. 


INTRODUCTION 


The purpose of this report is to describe and present instructions 
for using the computer programs developed by NASA/Ames Research Center 
personnel for calculation of the flow field about space-shuttle vehicle 
(SSV) configurations (refs. 1 to 3). The programs as described in this 
report are coded in FORTRAN for the CDC 7600 computer. The main programs 
and most subroutines are in a form resulting from translation from the 
special FORTRAN-like language (CFD) developed at NASA/Ames Research 
Center for use on the ILLIAC IV computer. The use of the CFD language 
results in a FORTRAN code which is more efficient when run on the CDC 7600 
than the pure FORTRAN code in which the programs were originally developed. 
However, the FORTRAN code as translated from the GDC language may prove 
difficult to read. A few subroutines such as input and output routines 
and a subroutine for calculating the body geometry have been coded directly 
in FORTRAN so that these programs cannot be run as they exist on the ILLIAC 
IV computer. The programs could be run, with minor modifications, on 
other serial processing computers such as the IBM 360 or other computers 
with sufficient memory capacity. The largest of the five programs requires 
approximately 650,000 octal words (220,000 decimal words) of core storage. 

A total of seven permanent tape or disc files are required by the programs 
for input and output in addition to punched card input data. 

The programs described in this report are • 

1. Finite-Volume Time Marching Code (FVTM) 

(ref. 1) 

2 . Finite— Volume Space Marching Code ^(FVSM) 

(ref, 2) 

3. Finite-Difference Space Marching Code (FDSM) 

(ref. 3) 

4. Streamline Coordinate Metric Coefficient Code (CHAOS) 

(ref. 4) 

5. Coordinate Transformation Interface Code (CTI) 
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The first program provides starting conditions for the second program 
which in turn provides starting conditions for the third program. The 
fifth program transforms the results of the first two programs into the 
coordinate system and computational mesh required by the third and fourth 
programs . 

The next sections of this report will describe in more detail the 
general use of the programs and then present instructions for preparing 
the input data and interpreting the output 

description of programs 

Because of the rather complicated nature of the flow about the space- 
shuttle vehicle, several individual codes have been developed rather than 
one overall code in order that the several different regions of the flow 
which occur can be handled with numerical efficiency. The three regions 
of applicability of the programs described in this report are shown in 
the sTcetch below. 



The individual programs listed previously will now be described briefly. 

Finite-Volume Time Marching Code 

This program is based on a second-order accurate, time-dependent, 
finite volume numerical algorithm that solves the unsteady flow equations 
for a perfect gas or a chemically reacting flow in integral conservation 
law form (ref 1) . The program is used to determine the structure of the 
flow field and the distribution of chemical species in the region of the 
blunt nose of the vehicle It carries the calculation through the sub- 
sonic portion of the field and continues to points :)ust beyond the sonic 
surface where it terminates. This code is applicable to a variety of 
blunt-body shapes and flows at supersonic Mach numbers up to 22 and 
angles of attack up to 40°. The computational mesh for the program is 
determined by a series of nested cones as shown in figures 1 and 2. The 
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surfaces defined by the cones shown in figure 1 are intercepted by a 
plane rotated about the cone axis in specific angular increments as shown 
in figure 2. Finally, a third set of surfaces which divide the region 
between the body and the shock into a specific number of layers produces 
a set of finite volume hexahedra to which a finite— difference computational 
scheme is applied to solve the flow equations in integral conservation 
law form. The flow field is calculated as the asymptotic solution of an 
unsteady flow started from arbitrary initial conditions. 

Finite-Volume Space Marching Code 

The second code is based on a second-order accurate, steady, hyper- 
bolic finite-volime algorithm that solves the flow equations in integral 
conservation law form for supersonic flow (ref. 2) . This code uses 
initial data determined by the previous program which are supplied on a 
curved surface ^ust downstream of the sonic surface (fig. 3) . The basic 
computational surface is circular conical as for the previous program, 
defined at each step of the calculations by the location of the apex of 
the cone, the cone vertex angle and the cone axis orientation with the 
body axis. In the present version of the program, the location of the 
apex of the cone is held fixed while the cone vertex angle, 0, increases 
monotonically toward 90° and the axis orientation angle, approaches 
zero so that at the final step the surface degenerates into a plane 
normal to the body axis 

Finite-Difference Space Marching Code 

The third program is based on a second-order accurate, hyperbolic, 
finite-difference algorithm that solves the flow equations for a perfect 
gas or a chemically reacting flow in differential conservation law form 
for supersonic flow (ref. 3) This program accepts initial data from 
the previous program on a plane normal to the body axis and then carries 
the calculation back over the remainder of the vehicle or until the Mach 
number becomes less than one in the marching direction (fig. 4) . The 
governing equations are solved between the body and the outer most shock 
wave which is treated as a sharp discontinuity. Secondary shocks which 
form between these boundaries are captured automatically. 

Streamline Coordinate Metric Coefficient Code 

The fourth program is based upon the partial differential equations 
relating the metric coefficients for the three-dimensional streamline 
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coordinates to body geometry and the three-dimensional inviscid solution 
Consequently, this program requires input data from all of the previous 
three codes. It determines the surface streamlines and then by using the 
axisymmetric analogy, estimates the heat transfer ho the vehicle (ref. 4) 
This program is not described in detail in this users manual. However, a 
list of the punched card input required is included in the Appendix. 

Coordinate Transformation Interface 

A fifth program is used to transform the data from the first two 
programs to the coordinate system and computational mesh required by the 
third and fourth programs. The program employs linear interpolations of 
the data of the first two programs to obtain data at user- specified axial 
stations and equispaced values of the meridional angle, and the radius 
r (fig. 4) . The data on the body surface are adjusted to satisfy the 
conditions of constant total energy and flow tangent to the surface. No 
attempt is made to maintain the surface entropy at a constant va lue for 
the perfect gas case. 


PROGRAM ORGANIZATION 

In this section, the general organization of the program will be 
described Specific information on data required for input and data 
developed for output will be described in subsequent sections. 

A general flow chart of the overall program organization is shown in 
figure 5. Each program requires some punched card input and some input 
data from disc or tape data files. Each program in turn produces new 
disc or tape data files, and printed output. 

Program 1 (FVTM) requires, initially, data from cards describing 
the free— stream conditions and the computational mesh. In addition, 
input are required describing the vehicle geometry. It is most expedient 
to store the geometrical input data on a separate disc file since the 
same data are required by program 2 and by the interface program. The 
specific construction of this file will be described in a subsequent 
section. Program 1 produces as output lists of the appropriate flow 
field quantities plus a disc or tape data file (data file 1, PlDATA} 
which has three possible uses . 
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1. To restart program 1 for continuation of an unfinished calculation. 

2. To start program 2 (FVSM) 

3. As input to the interface program (CTl) 

Program 2 requires some card input data describing certain print 
options, certain output options, and telling the program whether it is 
starting from program 1 or continuing an unfinished calculation. Program 2 
also requires the same geometry input file as did program 1. Program 2 
produces printed output and two data files. One file (data file 2, P2RS) 
is used to restart program 2 in case of premature termination. That file 
IS replaced every tenth step of the integration. The other file (data 
file 3, P2DATA) is used by the interface program both for producing a 
combined file from the results of programs 1 and 2 for the entire nose 
region and to produce a data set to start program 3. Data file 3 con- 
sists of the flow variables, chemical spepies concentrations, computational 
mesh coordinates and certain parameters for every tenth step of the inte- 
gration. The interface code transforms the accummulated data from pro- 
grams 1 and 2 to a body-axis oriented cylindrical coordinate system and 
interpolates the data to planes normal to the body axis at specified 
axial locations. The program then produces a data file (data file 4, 

CHAOS) to be used by program 4, the steamline metric program, and another 
(data file 5, BBDATA) for program 3, the steady shock-capture code Pro- 
gram 3 continues the calculation from where program 2 terminated, calcu- 
lating downstream over the remainder of the vehicle. Program 3 produces 
a data file (data file 6) to be used to restart in the case of a premature 
termination. Data file 6 is not produced in a form that is immediately 
capable of being used by program 4 To do so will require development of 
an interface program to convert the data to that required by program 4 as 
IS done by program 5 for the data of programs 1 and 2. 

INPUT TO THE PROGRAMS 
Tabular Form 

The input data required for calculating the flow field about a space- 
shuttle vehicle consists of several punched cards containing parameters 
describing the free-streara flow conditions, the computational mesh, and 
certain options that are available in the programs. A dictionary of the 
input data is presented in the next section. The data required to describe 
the body geometry are included since the data must be prepared on punched 
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cards and then loaded on a disk or tape file for use by the programs. 
Figure 6 shows the input variables as they are to be punched on the data 
cards for each program Further explanation of the proper preparation of 
input data is given in the section "SPECIAL INSTRUCT lOISIS FOR PREPARING 
INPUT DATA. " 


Dictionary of Input Variables 

The variables required for input on punched cards are defined in 
this section in the order in which they are required Additional details 
on the format of the punched data are given in figure 6. 


Procrram 1 { p VTM) 


DISKIN 

DISKOT 

NDEND 

NPRNT 


NPRT 


NEND 


logical variable indicating whether input is from disk 
file (T) or cards (F) 

logical variable indicating whether an output disk 
file IS to be written (T) or not (F) 

integer indicating the number of iteration steps to be 
computed in present calculation 

integer indicating whether final flow field results 
are to be printed 

0 - no 

1 - yes 

integer indicating whether initial computational mesh 
parameters are to be printed 

0 - no 

1 - yes 

integer indicating number of iteration steps to be 
calculated in first calculation 


IE 


integer, total number of latitudinal mesh points, 
maximimn 16 


JE integer, total number of shock layer mesh points, 

maximum 192/IE 

KE integer, total number of meridional mesh points, 

maximum 20 


JSHK 


6 


integer indicating location of shock, usually 
JSHK = JE-2 



NGAS 


EMACH 

PINF 

RINF 

GAMMA 

SDO 


C2,S3 

ALPHA 


integer indicating type of gas being calculated 
NGAS = 0 , nonequilibrium air 
= -1, perfect gas 

free-stream Mach number, M 

free-stream static pressure, dynes/cm^ 

free-stream density, gm/cm® 

free-steam ratio of specific heats, T 

initial estimate of shoch standoff distance, nose 
radii (fig. 7), see "SPECIAL INSTRUCTIONS FOR 
PREPARING INPUT DATA" 

parameters describing shock shape, see "SPECIAL 
INSTRUCTIONS FOR PREPARING INPUT DATA: 

angle of attack, degrees 


For information on how to determine the following input quantities for 
program 1, see "SPECIAL INSTRUCTIONS FOR PREPARING INPUT DATA" 


ZFOCNL 

THWL 


THLL 


RBX 


axial location of vertex of largest cone of computa- 
tional mesh, nose radii (fig, 7) 

angle between the free-stream velocity and the wind- 
ward intersection of the largest cone of the computa- 
tional mesh with the pitch plane, degrees (fig. 7) 

angle between the free-stream velocity and a line 
drawn from the axial location of the first computa- 
tional mesh cone vertex denoting the point of inter- 
section of the leeward limit of the computational 
mesh with the body, degrees (fig. 7) 

axial location of first mesh cone vertex, nose radii 
(fig. 7) 


Program 2 ( fvs M) 


NRSTRT integer indicating source of input data 

0 - data are obtained from the file produced by 

program 1 (FVUBB, data file 1) 

1 - data are obtained from the restart file produced 

by program 2 (data file 2) 
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NPFLOW, NPAREA 
NPXYZjNPFLUX 


DTHETP 


integers indicating whether output is to be printed 
for flow variables, mesh areas, mesh coordinates and 
fluxes, respectively 

0 - no output printed 

1 - output printed 

angular interval for output, degrees 


Program 3 (FDSM) 


GRFGS 


MODINjLTIN, 
MODOUT, LTOUT 


logical variable indicating whether CRT display is to 
be used. 

FALSE - no 
TRUE - yes 

Note* The present version of program 3 does not 

include the capability of employing the CRT 
display so the variable GRFCS should always 
be input as FALSE 


integers indicating input and output mode, source of 
input data and disposition of output data. MODIN, 
MODOUT, LTIN and LTOUT are used to control the 
following options * 


MODIN = 1, LTIN = 1 read from disk, logical unit 20 

MODIN = 1, LTIN = 2 read from disk, logical unit 21 

MODIN = 1, LTIN = 3 read binary version of file 

"BBDATA" 

MODOUT = 2, LTOUT = 1 write on disk, logical unit 20 
MODOUT = 2, LTOUT = 2 write on disk, logical unit 21 
MODIN = 3 read cards, unit 31, LTIN ignored 
MODOUT = 4 punch cards, unit 30, LTOUT ignored 
MODOUT = 5 write on disk, unit 22, LTOUT ignored, 
continueous storage of results 
MODIN = 6 LTIN = 0 read from last entry of con- 
tinuous storage file, unit 22 


MODIN =6, LTIN > 0 read from the first entry of 

the continuous storage file 
at which the integration step 
number is greater than or 


equal to LTIN, unit 22 
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NIT 

MODCFO 

MODBSO 

CRASHZ 


integer, upper limit on number of integration steps 
for current calculation 

axial increment at whicb complete flow field is 
to be printed, cm 

integer niomber of integration steps between output of 
body and shock data 

axial increment between output for storage file (Data 
file 6) 


NUMDtFM 

ZRMSH 

BETTA, PHIZRO, 
AR, BR 

ZFACT 

FCTROW 

NMPTS 

ZDMPM 

DCM 

NRPTS 

ZDMPR 

DCR 


integer number of mesh data cards to be inputted 

axial locations at which the computational mesh will 
be revised, cm. See "SPECIAL INSTRUCTIONS FOR 
PREPARING INPUT DATA- Clustering Parameters". 

parameters for clustering the computational mesh, see 
"SPECIAL INSTRUCTIONS FOR P3REPARING INPUT DATA - 
Clustering Parameters". 

axial locations for changing weighting factor on 
integration step size, cm, see "SPECIAL INSTRUCTIONS 
FOR PREPARING INPUT DATA - Step Size for Program 3" 

weighting factor on integration step size, may range 
from 0. 5 to 1.2 depending on rate of change of solu- 
tion 

integer number of data cards to be inputted for merid-t- 
lonal damping 

axial location at which meridional damping factor 
will be changed, cm 

meridional damping factor, see "SPECIAL INSTUCTIONS 
FOR PREPARING INPUT DATA - Damping Parameter" 

integer number of data cards to be inputted for 
radial damping 

axial location at which radial damping factor will be 
changed, cm 

radial damping factor 
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Program 4 (CHAOS } 


A dictionary of input data for program 4 taken from reference 4 is 
presented in the Appendix, 

Program 5 (CTI) 

JLIMIT integer, nvimber of shock layer locations at which 

data are to be calculated for program 4 (CHAOS) , 
including the body surface, but not the shock wave, 
maximum 24 

JL3 integer, number of shock layer locations at which 

data are to be calculated for program 3 (FDSM) , 
including the body surface and the shock wave, 
maximum 22. 

NUMK integer, number of meridional (i>yJ locations at which 

data are to be calculated for programs 3 and 4. May 
be different from (KE-1) used in programs 1 and 2^ 
maximum 20 

lOPT integer, indicates which output is to be produced 

3 - output IS produced only for program 3 

4 - output is produced only for program 4 

43 - output IS produced for both programs 3 and 4 

If any number besides one of these three is inputted, 
the program stops after printing an error message. 

NPRNT integer, indicates whether output of the calculated 

flow field data produced for programs 3 and 4 is to 
be printed 

0 - no 

1 - yes 

IZ integer, total nxunber of axial stations at which data 

are to be requested for program 4, maximum 40 

ZMTRC(I) array of axial stations at which data are requested 

for program 4, in fractions of the length (ZMAX) of 
the region covered by programs 1 and 2. Up to 40 
values are allowed. Values larger than 1. 0 and 
smaller than Z (1, 1, 1) /ZMAX of the computational mesh 
will be rejected. 


10 



Body Geometry Data 


"See "SPECIAL INSTRUCTIONS FOR PREPARING INPUT DATA - Body Geometry" 


RNOSE 
CONANG 
NSC (N) 

NS 

K 

NCF (NS,N) 

YPl 

YP2 

ZL(NS,N, 1) 

ZL(NS,N,2) 

CF (I,NS,N) 
1 = 1 , 2 , 

NCF (NS, N) 


radxus of nose, xnches 

angle of nose cone frustum, degrees 

integer, number of segments used to define curye fit 
for variable N, maximum NSC(N) = 13 

integer, segment number corresponding to a particular 
set of coefficients of a curve fit 

integer, order of the polynomial used for the curve 
fit, maximum K = 9 

integer, number of coefficients in curve fit for 
variable N in segment NS, maximum NCF(NS,N) = 10 

slope at forward end of segment NS 

slope at aft end of segment NS 

axial location of forward end of curve fit for 
variable N, in segment NS 

axial location of aft end of curve fit for variable 
N, in segment NS 

coefficients of curve fit for variable N in segment 
NS 


SPECIAL INSTRUCTIONS FOR PREPARING INPUT DATA 


Units 

Nondimensional quantities have been used for the input data where 
convenient. However, the free-stream density and pressure and the quan- 
tities describing the body geometry must be input in dimensional units. 

The user ' s attention is called to the fact that the FREE STREAM CONDITIONS 
are required in CGS UNITS while the GEOMETRY PARAMETERS are input in 
ENGLISH UNITS (inches) Angles are input in degrees. 

Shock Shape Parameters 

The initial estimate of the shock shape is specified by the three 
input parameters SDO, C2 and S3. The initial shock surface is generated 
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by a quadratic function of the latitudinal angle, 0, and the meridional 
angle, 0 (fig. 2) and is positioned at an estimated standoff distance, 
A^. The input parameter SDO is the nondimensional standoff distance 


SDO 



where is the nose radius. Experience has indicated that the initial 

estimate of the standoff distance should be larger than the converged 
value. Exactly how much larger it must be and what is the exact effect 
on the convergence has not been thoroughly investigated. For spheres, 
the shoch standoff distance is a weak function of Mach number for hyper- 
sonic Mach nitmbers. For bodies with spherical noses, a reasonable value 
for SDO IS 


SDO 0. 25 

Different values may be used if desired for specific cases 

Since the real shock is not expected to conform to the shape of the 
body, the remaining parameters are used to describe the distance between 
the shock and the body surface along rays of the computational mesh 
(fig. 2). Thus, the following expression describes the shock shape, 

S^k = SD0ri + 0f]^(C2) + 0^^(2 + cos 0^)S3] 

The parameters C2 and S3 can be adjusted to provide a desired degree 
of expansion of the shock layer. Values of C2 and S3 of 0.5 and 0.0 
respectively have been found to yield good results for bodies with sphere- 
cone noses Those values are used in the sample calculation to be pre- 
sented subsequently. Those values should also yield good results for other 
cases. Nevertheless, the quantities are required as input in order to 
provide the user with the ability to adjust the initial conditions. 

Computational Mesh Parameters 

The computational mesh is constructed from the input values of 
ALPHA, ZFOCNL, THWL, THLL and RBX. The definitions of these quantities 
are illustrated in figure 7 where a sketch of a typical sphere-cone body 
at an angle of attack is shown. The various parameters are necessary in 
order to define the computational mesh to account for the asymmetry of 
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the flow. The objective is to define the computational mesh to enclose 
all of the subsonic flow in the nose region plus supersonic flow with 
Mach numbers less than approximately 1.2 Thus, the limits of the mesh 
should be in supersonic flow slightly downstream of the sonic surface. 

For a discussion of the asymmetry of the sonic surface, see reference 5. 
For a sphere, or a sphere-cone at zero angle of attack, the value of RBX 
should be 1.0, ZFOCNL should be 1,0, and the angles THWL and THLL should 
be equal and of the order of 50° for hypersonic Mach numbers. The exact 
value used can vary with the free-stream Mach number. The reader is 
again referred to reference 5 for examples of the shape of the sonic 
surface on spheres. 

For sphere-cone bodies at nonzero angles of attack, the value of 
RBX should be 1 0, and ZFOCNL should be greater than 1.0, the exact 
amount being dependent on the angle of attack. The value of ZFOCNL is 
also somewhat dependent on the value chosen for THWL, The sum of THWL 
and ALPHA must not exceed 90° in order to allow program 2 (FVSBB to 
operate. A rule of thumb is to. limit the values so that 

THWL + ALPHA <80° 

This means that for angles of attack less than about 30°, THWL can be 
specified as 50° as for a sphere but should be decreased slightly for 
higher angles of attack. 

The values of ZFOCNL and THLL must be increased as angle of attack 
IS increased to accommodate the downstream shift of the sonic surface on 
the windward side of the body as described in reference 5. For THWL 

around 50° a first approximation to ZFOCNL can be obtained from figure 8. 

Since the exact shape of the sonic surface is dependent on the body shape 

and the free-stream Mach nimber as well as the angle of attack, the value 

obtained from figure 8 may require some adjustment for specific cases. 

f 

On the leeward side of the body, the sonic surface is usually situ- 
ated such that the parameter THLL can be specified in the range 50° to 
90°. 

A test IS made in program 1 at the end of each iteration to determine 
if any Mach numbers less than 1. 2 occur at the exit boundaries of the 
computational mesh. If any low Mach number flow is present, the user is 
so advised by a printed message "LOW MACH NUMBERS AT EXIT BOUNDARY. " If 
the condition persists for many iterations, the computational mesh should 
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be expanded by adjusting the parameters ZFOCNL, THLL, or THWL. The calcu- 
latxon must then be re- initialized with the quantity DISKIH = F and the 
remaining cards of figure 6 (a) supplied as input. The values of the 
meridional index, K, at which subsonic Mach numbers were found are listed 
along with the above message as well as the number of the particular 
layer in the shock layer in which the Mach numbers were found. These 
values aid in determining how the input parameters should be changed. 

If low Mach numbers are found for low values of K, the mesh needs expand- 
ing on the leeward side of the body. If the low Mach nimbers are found 
for high values of K (near K = KE) , then the mesh needs expanding on the 
windward side. In either case, developing the optimum mesh may require 
adjustment of all three of the parameters ZFOCNL, THWL, and THLL, 

Interface Input Data 

The interface program (CTI) is designed to produce disk or tape data 
files for use by programs 3 and 4 (fig. 5) . The interface program con- 
verts the data in the coordinate systems of programs 1 and 2 into the body 
axis oriented coordinate systems of programs 3 and 4 and then interpolates 
to calculate the flow properties at equispaced meridional angles and at 
equispaced radial positions across the shock layer. The number of merid- 
ional angles at which data are to be calculated, NUMK, is arbitrary, up 
to a maximum limit of 20. However, if data for programs 3 and 4 are being 
calculated in the same job, both use the same value of NUMK. On the other 
hand, the number of radial positions across the shock at which data are to 
be calculated may be different for each program. If only the surface 
streamline data are required for program 4, JLIMIT may be input as 1. 

If more values are desired, JLIMIT may be as large as 24. For program 3, 
the number of radial positions is specified by the parameter JL3 This 
parameter may be as large as 22. 

If data are being calculated for both programs 3 and 4 (lOPT = 43), 
or for program 4 alone (lOPT = 4), the final input for program CTI is the 
number of axial stations at which data are to be calculated and a list of 
the required stations. The required stations are input as fractions 
of the length of the nose region. The length of the nose region is calcu- 
lated in the program from the coordinates of the computational mesh. The 
final value of the list should be 1. 0 if data are to be calculated for 
program 3. 


14 



If data are being calculated for program 3 alone (lOPT =3), the 
input for the interface program consists entirely of the single card con- 
taining JLIMIT, JL3, NUMK, lOPT, and NPRNT. In that case, it is auto- 
matically assumed that the required data will be calculated at the axial 
station corresponding to the end of the nose region (ZMTRC = 1.0). 

Body Geometry 

The simulation of an actual aircraft or spacecraft geometry is done 
by analytical approximations. A typical cross section of the body at a 
given longitudinal coordinate z is composed of the analytic functions 
illustrated in figure 9. These six segments require specification of the 
eleven geometrical parameters listed in the figure as a function of the 
longitudinal distance z. To completely describe a typical vehicle, some 
segments can have zero length. 

The longitudinal variation of the eleven geometrical parameters and 
their derivatives with respect to z are required by the program. The 
analytic approximation of the body is arbitrarily required to be single 
valued, continuous, and have continuous first derivatives. To ensure this 
requirement, discreet values of each of the parameters and their first 
derivatives are obtained at various stations along the body from working 
drawings of a particular vehicle or spacecraft configuration (these para- 
meters need not be defined at the same z coordinates) . Then for each 
of these parameters, a cubic polynomial which is determined so as 
to maintain the continuity requirements at the end points, is used to 
describe the variation between the points. 

Note that it IS assioraed that the curve fits are obtained for the 
various length parameters in ENGLISH UNITS (inches) . 

Input data forms for the body geometry data are presented in fig- 
ure 6(e). The first card required contains the nose radius and the nose 
cone angle (card no. 1 in fig. 6(e)). That card is followed by eleven 
groups of cards representing the curve fits for the eleven variables 
shown in figure 9. The first card in each group contains the variable 
NSC, the number of segments in the curve fit for that variable (card no. 2 
in fig. 6(e)). That card is followed by two or three cards for each seg- 
ment, The first card of each segment contains the information defining 
the curve fit in that segment, the slopes of the variable at the ends of 
the segment, and the axial locations of the ends of the segment (card 
no, 3 in fig, 6(e)). The second card contains the first five coefficients 
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of the curve fit in the segment (card no. 4 in fig. 6(e)). The third 
card contains the remaining coefficients if a high order polynomial is 
being used. Only four coefficients are used in the sample case discussed 
herein so only one coefficient card is required for the sample case for 
each segment. Up to thirteen segments may be used to described the varia- 
tion of a particular variable over the total length of the vehicle. A 
sample input geometry card deck is listed in Table I. A computer program 
for obtaining these data is described in reference 6. 


Clustering Parameters 

In program 3 (FDSM) the radial coordinate is normalized to allow for 
a reclustering of computational mesh points in the meridional direction 
at specified axial stations, ZRMSH (see Dictionary of Input Variables, 
Program 3) The reclustering is accomplished by introducing the inde- 
pendent variable transformations given by (ref. 3) 


z 


z 


i iz, T,<l>) = (r - r^)/{r^ - r^) 


H(0) = irjB + sinh“^,[(^^ - l^sin(B)J | /p 
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1 + (e*^ - 1)0 /tt 
o 
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n(0) = 0 


p = 0 


where p and 0^ (BETTA and PHIZRO) are arbitrary parameters that control 
the degree and location of clustering. No clustering occurs when p is 
zero. As p increases, the degree of clustering increases with the 
greatest density of points appearing about the ray defined by the angle 



Another set of parameters controls the reclustering of the computa- 
tional mesh points in the radial direction. The equations for this trans- 
formation are (ref 7) 
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I = a + (1 - a) 

Thxs transformation permits the mesh to be refined near the body alone 
(a = AR = 0) or to be refined equally near both the body^and the shock 
(a = AR = 0.5), The effects of different values of a and P (AR and BR) 
are shown in figure 10. 

Step Size for Program 3 

In program 3, the length of the step to be taken at each integration 
is determined to satisfy the mathematical requirements for stability. 
However, experience has shown that the mathematical requirements alone 
are not sufficient for accurate and efficient operation of the program. 
Therefore, a weighting factor (FCTROW) is provided so that the user may 
vary the step size. In regions where steep gradients are to be expected 
such as around canopies or wings, or other regions where the body radius 
may change rapidly with axial distance, the step size should be even 
smaller than dictated by the mathematical stability criteria (FCTROW < 1.0). 
On the other hand, in regions where the solution changes gradually the 
step size can be even larger than the stability criterion allows 
(FCTROW > 1.0) The exact variation of the weighting factor (FCTROW) that 
can be allowed must be determined by experience. 

Damping Parameters 

In order to avoid instabilities due to severe ‘pressure changes in the 
numerical solution of the finite-difference equations, program 3 employs 
a fourth-order damping terra which requires input of the axial location 
where the damping is to be changed and the new value of the damping coef- 
ficient to be used. The use of such damping is described in reference 8. 
Damping coefficients are input for the meridional and the radial direct- 
tions. The exact values to be used must be determined from experience. 
According to reference 8, destabilization is avoided if the damping 
factor IS less than or equal to 0,5. 

Restart 

Due to the long times required, programs 1, 2, and 3, have the 
capability of restarting at an intermediate step in the calculation. 
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The general procedure for running the programs is the same for restarting 
as for the initial start. However, some of the input data may require 
changes if the calculation is being restarted. 

In program 1, only the first card shown in figure 6 (a) is needed for 
restarting, but the disk or tape file created by the previous unfinished 
calculation must be made available for input. A file of the required data 
IS created every 25 iteration steps. The input quantity that must be 
different from the initial calculation is the logical variable DISKIN 
which is input as F for input of the initial data from cards and T 
for input of the restart file. In addition, the quantity NDEND must be 
input as the number of iteration steps to be calculated in the present 
calculation. (Note that NDEND is NOT the total number of steps to be 
calculated. The total number of steps is the number of steps that were 
calculated in the previous calculation, NEND, plus NDEND) . 

In program 2, both data cards shown in figure 6 (b) are always 
required. For starting initially, the first quantity on card no. 1, 

NRSTRT IS 0. For restarting, it is 1. The remaining quantities may be 
left the same or changed as desired. If DTHETP is changed on restarting, 
the printed output will be produced according to the previous value of 
DTHETP the first time and then will be produced at the new interval. A 
file of data for restarting is produced every ten integration steps. 

In program 3, only the first card shown in figure 6(c) is required 
for restarting. The quantity MODIN should be six for restarting. The 
quantity LTIN can be arbitrarily specified. As indicated in "Dictionary 
of Input Variables V if LTIN is zero, the input data will be read from the 
last entry of 'the storage file from logical unit niomber 22. Otherwise, 
if LTIN IS greater than zero, the data will be read from the first entry 
of the storage file at which the integration step number is greater than 
or equal to LTIN. 


OPERATING PROCEDURE 

In this section, the construction of card decks for operation of the 
computer programs is described. First, a general description of the 
operations required is given. Then the specific Job Control cards needed 
for operation on the CDC 7600 computer at NASA/Ames Research Center are 
listed. 
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General Job Control Sequence 

The following list is the general Job Control procedure that would 
be required to run programs 1, 2, 3, and 5 in sequence without stopping or 
restarting. The reader is referred to figure 5. 

1. Attach geometry input file 

2. Load program 1 (FVTM) 

3. Execute program 1, providing required card input data 

4. Store data file 1 produced by execution of program 1 

5. Load program 2 (PVSM) 

6. Execute program 2, providing required card input data 

7. Store data file 2 produced by execution of program 2 

8. Store data file 3 produced by execution of program 2 

9. Load program 5 (CTI) 

10. Execute program 5, providing required card input data 

11. Store data file 4 produced by execution of program 5 

12. Store data file 5 produced by execution of program 5 

13. Load program 3 (FDSM) 

14. Execute program 3, providing required card input data 

15. Store data file 6 produced by execution of program 3 

For continuing an unfinished calculation, the Job Control Sequence 
would be basically the same as in the previous list except that the 
appropriate restart data file would be attached and the part of the 
sequence that was already completed would be omitted with the exception 
that the geometry input file would be attached as before (step 1) . 

Logical unit numbers required for input and output are given names 
to avoid confusion in setting up the Job Control Card deck These names 
denote the files actually used by the programs and are listed in Table II. 
Permanent storage files may be labeled at the user's convenience. Normally, 
tape 5 is used for input from ptuichdd cards and tape 6 is used for printed 
output. In program 3, tape 31 may be used for punched card input and 
tape 20 IS used for punched card output. Other logical unit numbers are 
used for input or output from the permanent disk or tape files. 

Job Control on the CDC 7600 

Sets of Job Control Cards required to use the computer programs on 
the CDC 7600 computer at Ames Research Center are listed in Tables III, 

IV, V, VI and VII. It is assumed that the programs in binary form will 
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be stored on a disk or tape file. The reader is referred to the CDC 
SCOPE 2 Reference Manual for explanation of specific parameters used on 
the cards. 

In Table III, the Job Control cards necessary to run program 1 from 
an initial start are listed. The usual procedure for running program 1 
will be to do a complete calculation in at least two steps. The first 
step might run only a few iterations to verify that input quantities are 
correct and that the solution may be expected to converge. The next 
step would restart from the previous step and run a larger number of 
iterations. Criteria for convergence of the solution are discussed in 
the next section. 

In Table IV, the Job Control cards necessary to restart program 1 
and then run all programs except program 4 in a single 30 b are listed. 
Such a 30 b could require more than 1 hour for a perfect gas case, depend- 
ing upon how many iterations were required from program 1. Substantially 
longer times will be required if a real gas is calculated. Also, the 
storage file for program 3 (CONTOU) may become very large. If it is 
desired to store that file on a tape instead of disk storage, the cards 
32 and 33 should be replaced with the appropriate cards for transferring 
the file to tape. 

In Table V, the Job Control cards necessary to restart program 2 
and then run all remaining programs except program 4 are listed. See the 
previous discussion of Table IV for an alternate procedure for program 3. 

In Table VI, the Job Control cards necessary to restart program 3 
are listed. See the previous discussion of Table IV for an alternate 
procedure for program 3. 

In Table VII, the Job Control cards necessary to run program 4 
(CHAOS) are listed. 


Convergence Criteria 

In program 1, the calculation is controlled by the user-specified 
quantities NEND and NDEND, defining the number of iterations to be 
executed. No convergence criteria are built into the program. As an 
aid in determining the degree of convergence of the solution, several 
quantities are printed at each iteration. The most informative of those 
quantities are the x-coordinates and the pressure at the shock at points 
near the downstream boundary on the windward and leeward sides of the 
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body in th.e pitch plane. The coordinates appear as the quantities 
X(l,ILjJSHK) and X(KL, IL, JSHK) while the pressures are P(1,M), and 
P(KL-1,M) where M= (IL-2 )XJIh-JSHK- 1. (See the output described in the 
section entitled "NUMERICAL EXAMPLE"}. For a converging solution, those 
coordinates and pressures will reach a steady-state, nearly constant, 
value after many iterations. For a solution that is not converging, one 
or all of the values may oscillate v;ith significant, possibly increasing 
amplitude. Another quantity printed each iteration is the shock stand- 
off distance (printed as the z-coordinate of the computational mesh at 
K=l, 1=1, J=JSHK, Z (1,1, JSHK). That quantity alone is not a valid indi- 
cation of convergence since it converges* rapidly to a steady-state value 
while quantities in other parts of the computational mesh converge more 
slowly. Even more dangerous, it may appear to converge and then at the 
last step diverge suddenly. The shape of the shock near the downstream 
boundaries of the computational mesh appears to be a more sensitive indi- 
cation of convergence than the shape in the stagnation region. 

Since the total energy is not constrained in the calculation to be 
constant, the approach of the total energy to the free-stream value is 
another, perhaps the most important, indicator of convergence. A series 
of quantities related to the relative error between the calculated total 
energy at each point in the computational mesh and the freestream value is 
printed each iteration. The quantities are the maximum relative error, 
and the total numbers of errors less than 10 percent, 1 percent, and 0.1 
percent of the freestream energy. As for the shock standoff distance, the 
maximum energy error alone is not a sufficient indicator of convergence. 
However, a converging solution should have all the energy errors less than 
10 percent, most errors less than 1 percent and the number less than 0.1 
percent should show an increasing trend. 

At the end of a specified n\amber of iterations, the solution for 
all flow field quantities is printed. Included in that output is the 
Mach number distribution. That distribution is an excellent indicator 
of the quality of the solution since the Mach mmnber should increase 
monotonically and smoothly from the stagnation region (1=1) to the down- 
stream boundary (I=IL-1) in the latitudinal direction. 

Clearly, the determination of convergence is a subjective problem. 

The primary concern is that the flow field quantities and shock shape 
approach a steady state. Approximately 400 to 600 iterations are usually 
sufficienu to achieve a reasonable degree of convergence if the solution 
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IS converging, as indicated toy the total energy error. Examples of con- 
verging and nonconverging solutions are presented subsequently in the 
section "NUMERICAL EXAMPLE." 

MESSAGES PRINTED BY THE PROGRAM 

This section lists the messages given toy the programs and tells 
what to do when they are encountered. 

Messages from Program 1 

(1) READ FROM TAPE (11) ITERATION NO = 

This message is printed if program 1 is being restarted from a data 
set. It shows the iteration number at which the calculation is started. 

(2) WRITE TO TAPE (10) ITERATION NO = 

This message is printed every 25 iterations of program 1 as a 
restart data set is written to the disk (or tape) file associated with 
the unit ntmtoer 10, 

(3) NERR=n IN SERCH 

This statement may toe printed toy the subroutine GE0M3 during a 
calculation of the body surface gwometry if an error condition is encoun- 
tered. It probably means that one or more of the curve fits representing 
the body geometry was input in error. Check the body geometry input data 
or the quantities ZFOCNL, THWL, THLL, or RBX. 

(4) WAVE ANGLE LESS TEIAN MACH ANGLE 

This statement is printed if the shock wave angle is less than the 

Mach angle. It is most likely to occur at the initial step for low Mach 

numbers. Check the geometry input data and the input shock shape para- 
meters . 

(5) THIS WAVE IS A MACH WAVE 

This statement is printed if the shock angle is equal to the Mach 
angle. Check the geometry input data and the input shock shape para- 
meters, 

(6) ^ ERROR IN SHOCK 

This statement is printed if difficulty is encountered in calculating 

the pressure at the shock wave at the initial step. Check the geometry 

input data and the input shock shape parameters. 
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(7) MORE THAN 200 ITERATIONS 

This message is printed if the initial condition subroutine START in 
program 1 is unable to locate a computational mesh point compatible with 
the input geometry parameter curve fits . Check the geometry input and 
the computational mesh input. Note carefully the units required for 
dimensional quantities. 

(8) ***LOW MACH NUMBERS AT EXIT BOUNDARIES*** 

J=n K=n , n , n ... 

This message is printed if the computational mesh does not enclose 
all of the flow with Mach numbers less than 1.2, The points in the last 
cone of the computational mesh at layer J=n are listed to aid in deter- 
mining how the mesh should be modified. If low values of K are printed, 
the mesh needs expanding on the leeward side of the body. If high values 
of K (near K=KE) are printed, the mesh should be expanded on the wind- 
ward side. In either case, the best mesh may require adjustment of all 
three of the parameters ZFOCNL, THWL, and THLL. 

Messages from Program 2 

(1) ***ERROR- INTEGRATION STEP ZERO OR NEGATIVE 

This message is printed if the subroutine which calculates the 
integration step size produces a negative or zero value. One possible 
cause of such a condition is the occurrence of unity or subsonic Mach 
numbers in the calculated flow field. Such a situation can occur if the 
computational mesh of program 1 did not enclose all the subsonic flow. 

The situation can also occur when the solution obtained from program 1 
IS not adequately converged or does not extend sufficiently far into the 
supersonic region. In that case, program 2 may integrate several steps 
before the error condition is encountered. The probable remedy in all 
these cases is to expand the computational mesh and rerun program 1. 

Messages from Program 3 

(1) DATA READ FROM CONTINUOUS FILE-UNIT 22? NSTEP=n 

Z = . . , 

This statement is printed when the data have been read from the 
restart file at step n. This corresponds to the input option MODIN=6, 
LTIN=n. 
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(2) data written on continuous file-unit 23? 

NSTEP=n, Z= 

This statement is written each time the data are written to the 
storage file (every CRASHZ cm, see "Dictionary of Input Variables”). 

Messages from Program 5 

(1) INPUT ERROR lOPT MUST BE 3 OR 4 OR 43 

See Dictionary of Input Variables for Program 5. 

(2) UNSTEADY BLUNT- BODY TAPE READ, n STATIONS ON IT. 

This statement is printed when the file from program 1 has been 
completely read. The number of stations refers to the nimnber of computa- 
tional mesh cone surfaces. 

(3) STEADY BLUNT-BODY TAPE READ, n STATIONS ON IT. 

TOTAL NO OF STATIONS=m, 

This statement is printed when the file from program 2 has been 
completely read. The number of stations refers to the number of times 
data were stored in the calculations done by program 2. The total is 
n plus the previous number of stations read from program 1. 

(4) SUBSCRIPT LIMIT MISMATCH. KEMl=k, KLMl=n, I2=m 

This statement is printed after comparing the value of KEMl from 
the program 2 interface file (data file 3) with that of KLMl from data 
file 1. These values are the number of mesh cells in the meridional 
direction. The two values should be the same. If not, the two data 
sets being read are not compatible. The value of 12 is the total number 
of computational mesh points at each meridional position. 

(5) SUBSCRIPT LIMIT MISMATCH JSHK2=k, JSHK=n, I2=m. 

This message is printed if the number of computational mesh points 
across the shock layer do not agree between data sets 1 and 3 (fig, 5) . 
See message (4) . 

(6) I2=n. TOO BIG 

This message is printed if the total number of mesh points in the 
shock layer at each meridional position for the combined data sets 
produced by the data from programs 1 and 2 exceeds 300. The remedy is 
to redo the calculation, using a coarser grid of points in program 1 by 
either using smaller values of JSHK and/or IL or by expanding the 
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computational mesh, so that the last mesh cone is mote nearly planar and. 
noimal to the body axis. 

(7) MORE THAN 40 AXIAL STATIONS REQUESTED EXCESS IGNORED 

This statement is printed if there are more than 40 values in the 
input array of axial stations, ZMTRC(l), at which data are to be calcu- 
lated for program 4 

(8) ALL Z-LOCATIONS. GT. NO. n IGNORED. TOO BIG. 

This statement is printed if some of the input z values are aft 
of the most aft station of the computational mesh 

(9) THE FIRST n Z-LOCATIONS IGNORED AS TOO SMALL 

This statement is printed for those input z stations which are 
forward of z of the computational mesh. 

(10) ERROR IN BODY SURFACE DATA STATION . . . DISREGARDED. 

This statement is printed whenever the interpolation and extra- 
polation of data from programs 1 and 2 is unable to satisfy the condi- 
tions of tangency and constant total energy. That problem may arise at 
stations near the nose where certain points of the computational mesh 
would have small velocities which must be calculated as a difference 
between large approximate quantities. In such a case, the calculation 
may produce a negative result from the expression 

p) 

when such an error occurs the data for the z station is not printed 
and IS not placed on the data set. 

NUMERICAL EXAMPLE 
Input Data 

A list of the punched card input data for a sample calculation on 
a space shuttle configuration is shown in figure 11(a), The case being 
calculated is for a Mach number of 10, angle of attack of 30° at an 
altitude of 50 kilometers in air treated as a perfect gas with ratio of 
specific heats equal to 14. 

The punched card geometry input data are listed in Table I. Those 
data are stored on a disk or tape file as described previously. 
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Printed Output 


As can be noted in the input shown in figure 11, all output option 
parameters were input so that all output would be printed. Only selected 
pages of the output will be presented herein. Figures are presented 
subsequently to illustrate the results of the calculations. 

Procrram 3 - The first printed output (fig. 12) from program 1 is a 
list of the important input quantities plus other free-stream quantities 
calculated from the input data Those additional quantities are* , 


VINE 

free-stream velocity 

AINF 

free stream speed of sound 

HINF 

total energy 

RMZ 

P V 

EINF 

total energy per unit volume 

EIINF 

specific internal energy 

TINF 

free-stream temperature 


The next list of output is a table of the x-coordinate, RFOC(K, 1,1), 
and the z-coordinate, RFOC(K, 2,1), of the computational mesh cone vertices 
listed with the meridional variation, K, down the page and the axial 
variation, I, across the page. Next, a table of the body radii associated 
with the rays of the computational mesh is printed with the meridional 
variation, K, down the page and the axial variation, I, across. These 
tables may be omitted by letting NPRT=0 

The next output is printed at the first integration step. It is 
basically a reprint of the input data with dimensionless input quantities 
converted to dimensional quantities. 

The next output is always printed. No options are provided for 
omitting it Two groups of data are printed at each iteration step to 
aid in monitoring the progress of the iteration. 

The first group of data are values of eigenvalues and other para- 
meters related to the numerical solution of the finite— difference equa- 
tions: The x-coordinates of certain points in the computational mesh 

and the pressure at those points are printed along with the shock stand- 
off distance, Z(1,1,JSHK). These quantities have been discussed previously 
in the section entitled "Convergence Criteria. " 
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The second group of data are a measure of the accuracy of calculation 
of the total energy. The maximum total energy error relative to the free- 
streara total enthalpy is printed along with the indices, i, 3 ,k, of the 
computational mesh at which it occurs. In addition, all the energy errors 
are counted and the number of errors less than 10 percent, 1 percent, and 
0.1 percent are printed. The total number of errors being counted is 
( (IL-1)X(JL-1)X{KL-1) ) . 

The remaining output from the program is selected with the input 
parameter NPRNT on the first input data card. First is a list of all the 
energy errors relative to the free-stream total energy. They are listed 
in groups of (IL-I)X (JSHK-1) values for each meridional position, K. 

Next IS printed the flow variables, p, u, v, w, p, Mach number, and 
e and, if a chemically reacting gas is being calculated, the chemical 
species concentrations. Finally, the coordinates of all the computational 
mesh points are printed. 

If NPRT and NPRNT are input as zero, only the input data and the 
Iteration step data are printed. 

The results produced by program 1 for the sample case are presented 
graphically in figures 13 through 17. In figure 13, the shock stand-off 
distance is seen to converge rapidly to a steady value within 200 itera- 
tions while the maximum energy error (fig 14) takes much longer. The 
energy error count shown in figure 15 indicates a converging solution with 
the number of small errors increasing with iteration. The shock shape 
shown in figure 16 is smooth at 200, 400 and 600 iterations, but moves 
slightly as the calculation proceeds. Clearly in this particular case, 
the shock shape alone is not a sufficient indicator of convergence since 
It is smooth and contains no perceptible perturbations even after only 
200 Iterations. A better indicator of the quality of the solution is the 
Mach number distribution, a sample of which is presented in figure 17. 

In that figure, the Mach number distribution at the shock in the plane 
of symmetry is plotted as a function of the nondimensional latitudinal 
index, I. At 200 iterations considerable roughness exists in the distrib- 
utions which IS gradually smoothed out as more iterations are calculated. 
Finally, the output quantities which seem to best serve to indicate the 
trend of the solution are shown in figure IS. Those quantities are the 
x-coordinates of the ends of the shock in the symmetry plane. The shock 
IS seen to initially oscillate markedly on both the windward and the lee- 
ward side but to eventually reach a steady state. 
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To illustrate the behavior of the calculated results when convergence 
IS not achieved, figures 19-24 are presented. The calculations are for 
the same sample case as figures 13-18, but with ZFOCNL=5.0. The most 
notable features of these figures are that the solution for early itera- 
tions follows similar trends to those of the convergent solution shown 
previously. In particular, the shock stand-off distance (fig. 19) con- 
verges rapidly in both cases and remains converged. Other quantities 
appear to converge and then diverge. 

Procfram 2 - The first output from program 2 (fig. 25) is the itera- 
tion step of program 1 at which the solution was obtained followed by a 
list of the parameters and constants obtained from program 1. Included 
in that list are the location of the vertex of the computational cone 
(x^,z^), the initial angle of the cone axis from the horizontal (z) axis, 
BETAO and the initial cone half angle, THETAO. Other data in the list 
are the same as output initially from program 1. 

The initial output is followed, if NPFL0W=1, by the list of densities, 
velocities, pressures and Mach mmibers obtained from program 1. However, 
the quantities in this list do not correspond exactly to values from 
program 1. A linear interpolation is performed to go from the finite 
volumes of the program 1 computational mesh wherein the flow quantities 
are known at the center of the mesh hexahedra to the finite areas of the 
program 2 computational mesh wherein the quantities are known at the 
center of the finite areas. 

Since the total energy is not constant in the solution from pro- 
gram 1, an adjustment is made in the initial data of program 2 to ensure 
that the total energy equals the free-stream value throughout the flow. 

The next output is the adjusted initial data, followed, if NPFLUX=1, by 
the calculated initial flvixes of mass and momentum. After the fluxes 
are calculated, the flow quantities are calculated and printed to verify 
the accuracy of the decoding procedure. 

The next output, if NPXYZ=1, is the list of mesh coordinates at the 
step, followed immediately by those at the (N + 1)'’^^ step. Next, 
if NPAREA=1, IS a list of the areas of the projections of the mesh cell 
faces on the coordinate planes at the initial step. 

Subsequent output is dependent upon the quantity DTHETP as well as 
the option parameters. In the case shown in figure 25, the output was 
required every 5°. At every integration step, the current value of the 


28 



cone angle, THETA and the integration step size are printed along with a 
list of the body and shock radii at each meridional position. Every 
DTHETP degrees , other data are printed as determined by the input option 
parameters . 

Program 5 - The interface program prints first the free-stream 
quantities, and the input and adjusted arrays of the axial stations at 
which data are to be calculated from programs 3 and 4. In the example 
shown, the input stations were all within acceptable limits, so the adjusted 
list of axial stations is not printed. Next, the values of the flow 
quantities which are calculated at each z station (fig. 26) are printed. 

In addition, a list of the body radii as computed from the exact relations 
and from interpolation is printed at each station, along with the relative 
error, the derivatives of the body radius with respect to z and 0 and 
the shock radius. When all the data for program 4 have been printed, a 
message is printed, saying 

***CHAOS TAPE WRITTEN*** 

This refers to the generation of data set number 4 in figure 5. Next, 
the data for progrcim 3 are printed. Note that the output for program 4 
IS in ENGLISH UNITS (inches) while that for program 3 is in CGS UNITS as 
for programs 1 and 2. The required unit conversions are done in the 
interface program (CTI) . Since program 3 makes use of the shock radius 
and its derivative with respect to z and <p, those quantities are added 
to the final output list' The final output is a message saying 

***SHOCK-CAPTURE CODE STARTING TAPE WRITTEN*** 

which refers to the generation of data set number 5 in figure 5. That 
data set contains the results of the blunt-body calculations to initial- 
ize the calculation for the remainder of the vehicle. 

Program 3 - The first page of output from program 3 (fig. 27) is a 
list of the quantities input from punched cards. This is followed by 
a list of some of the quantities provided by the data set produced by 
program 5. Next, the calculated quantities for the entire computational 
mesh at the initial axial station are printed. After that, the flow 
quantities are printed for the body surface and the outer shock at every 
second integration step. 
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Additional Calculations 


In addition to the sample case discussed herein, perfect gas calcu- 
lations have also been made with programs 1, 2, and 5, for Mach numbers 
of 5.0 and 7.0 at an angle of attack of 30° and for Mach numbers of 5.0, 
7.0, and 10.0 at an angle of attack of 35°. A calculation for a real gas 
at a Mach number of 27.5 and an angle of attack of 41.4° has also been 
completed. The results appear to converge in all those cases in the same 
manner as described herein. The upper and lower limits of applicability 
of the programs with regard to Mach number and angle of attack have not 
been investigated. 
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TABLE I 


BODY GEOMETRY INPUT DATA 


28 . 

- 11 

I ^ -C.2144500E 01 -0.94895O9E 00 0.2620000E 01 0.7250000E 01 

-n.i-iqiOOOE ''2 -0.214A500F 01 0.1978716E 00 -0 . 990 1 192 E-02 

7 3—' a -0.94S9599E' 00- -0,"35A1200E“‘00 -0T725000OE 01 0.3500000E C2 

-0.18‘=OOOOF -0.948959°F 00 0 . 197559 lE-01 -0 . 2 1 80920 E-03 

3^4 -G.^'^41 200>^ 00 -0.2262800E 00 C.3500000E 02 0.6500000E 02 

-0.’4?5000F 02 -0.-3F4i200B~00 0-.-28 1 7^2-7E-=&2 -0,“15Y591'2E^04 

434 -0.2262800E 00 -0.1227600E 00 0.6500000E C2 0.1150000E 03 

-0.4275000E 0? -0.226?8nOF 00 0. 1606797E-02 -0.7623955E-05 

5-^4 -0«1227800E" 00 -0 . 7869995E---0-1 OrM-50000E 03 -0«1650000E 03 

-'^,3inn000'=^ ^2 -0.1227800E 00 0 . A851967F-03 -0. 5919579E-06 

6 4 -0.7»6a995E-01 -0 . 2000000E-0 1 0.1650000E 03 0.2650000E 03 

,'i6'^0''00P ^2 -*'^.78649<>5='-01 -OTl239992f’=^ Ovlrr3t)003E-O5 

7^4 -0,?^000noE-9l -0. 20000noE-01 0.265'^OOOE 03 n.6650000E 03 

-'^,8190000'=' 02 -0.?00''000F-C1 -0. 1 788139E-1C 0.298C232E-13 

a - -7 4 ^noaOn^F^Ol -0. 200C000E-01 0 ^ 665000 nF 03 'OlOOSOOCOE 03 

-n.AOS'^OOOE 02 -0.2000000E-01 -0 . 31 7891 1 E-1 0 0 . 7064249 E-1 3 

974 -0.2000000F-01 C.1256055E-01 0.9650000E 03 0.1072000E 04 

-0,7580009=^ ■''p-_o,- 2000000 F --01 0 , 5991 O 64 E-='- 07 , 0 r5'746'570E-06 

19 3 4 0. 1256055E-01 0 . 612COOOE-01 0.1072000E 04 0.1165000E 04 

~n.762'^C00F 02 0.1256055E-01 0 . 2563422E-04 0 , 1690452E-05 

n- ^ 4 0.6120000E--01 0 . 1 L 8 >^COOE- 00 0‘rl 1-65000E "04 0-rl2-93300E 04 

-O.7^FnO00‘^ 0? ' 0 . 6 I 2 OOOOF-OI Q.4043977E-03 -«0 . 95 1 1 1 32 E-06 — 

7 

1 4 ^,0000900- f'.-oaoofpo- •0.-Y6YOOOOE”01 0-,1-65COO'OE 03 

fi,ri''nnnon n.O'^onoon 0^00000 0*0000000 

274 9.0000000 0, 0001 COO 0.1650000E 03 0.1293300E 04 

r'.O'^O^OOO ^.0000009 ^,'^000000* '7-ooo'OOOQ 

11 

1 3 4 O.OOOOOOO O.OOOOOOO 0.2620000E 01 0.7250000E 01 

^,0009000 ^.00-0^6f>0- -0.01000000 OVOODOOOO 

234 0.0000000 -0.1033400E 00 0.7250000E 01 0.3500000E 02 

0.0000000 0.0000000 -0.2119718E-02 0.6191885E-05 

3 3- 4--0.1033400E- 00'--0.~l-2987-70E“00- 0ii350C00t)F 02 0.6500000E 02 

-O.IFroOOO'^ 91 -0.1033400E DO -0.4481035E-03 0. 1293452E-06 

434 -0.1298770E 00 -0 . 9276694E-01 0.6500000E 02 Q.1150000E 03 

-f',3nf)r>9nOF-91--i^rl-2'98Y-70F Tn- -nrr4958?8-P::-0-3 0:T694'243TE-05 

F 3 4 -0.9276694F-01 -0. 5375000E-01 0.1150000E 03 0.1650000E 03 

-O.nOOOnOF 02 -C.92756O4E-01 -0. 7443213E-03 0.1539321E-04 

4'-0.5-775UOOE'-01'--Or5375000E-0^1 OvlS-SOOOOE 03 0.66"50000E 03 

-0.1562800E 02 -0.5375000E-01 -0. 1831054E-09 0.2441406E-12 

774 -0.5375000E-01 -0. 1745500E-01 0.6650000E 03 0.7650000E 03 

-0.425000GE-0 2 — 0 .-5-3 7 50-0 OE-0 1 6 9 5-4 9i-E*- O 3 - 20 49 4 0 5 

834 -0.1745500E-01 0.4S03300E-01 0.7650000E 03 0.9650000E 03 

~O.450Ono0F 0 ? -0.1745500E-01 0. 3093849E-03 - 0 .485 5502E-06 

934 0.4807-300E^01 OrOOOOOUO 0V90-50 000 E"03- -OV107ZOOOE“04 
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TABLE I (Contmued) 


o,40o'ooocr 

02 

1 0 3 

4 

0 .•a4C'''000'=’ 

02 

11 ^ - 

4 

0 ,74omoo=’ 

o? 

12 


1 “3 

4- 

0,I187000‘=‘ 

02 

2 3 

4 

0 .1825000F 

02 

3 3 

4 

0 .7900000E 

02 

4 ^ 

4 

0,F2'iOOOOF 

02 

5 ^ 

4 

0,69«:rv0f'0- 

f>7 

6 

4 

0 ,8?'?'^000E 

02 

-t rt. 

’ 4 

0 .9500000E 

02 

S 3 

4 

0 .1 740000«^ 

03 

9 3 

4 

0 ,19O‘^000^ 

03 

to ' 3 

4 

n,9640000F 

03 

t ] 3 

4 

n ,4740000’=' 

0 7 

17 7 

4 

0 ,45«0000F 

03 

11 


1 3 

4 

O.ooooooo 


7 " 7 

4 

0.0000000 


7 7 

4 

0 .I'iOOOOOE" 

^1- 

4 7 

4 

0 .4OO0000F 

01 

5' 3 ’ 

4 

O.t IQOOOOP 

02 

6 3 

4 

0,~1 F^ptHlOP- 

02'" 

7 7 

4 

0,4250000'^ 

02 

8 "3 

4 

0.4300000F 

02 


0.4^0^30DE'-01 

0,0000000 

O.OOQOOOO 

0.000000'^ 

O.OOOCCOO 

0,-?T44500E 01 
0.2144500E 01 
O.IOOOOOOF 01 

o.iooooooE or 

0.5095200E 00 
O.SOO'SZOOE 00 
0-.39834‘90F 0-0 
0,^98^4<?°E 00 
0.291470QE 00 
0,?<514700E 00 
0.2171200F 00 
0.2171200E OO 
O,re4l00-0E'~OO 
0.1641000E 00 
0.1644000E 00 
0.1644000E 00 
0.4040000E 00 
0.4040000F 00 
O.IOOOOOOE OT 
O.IOOOOO^F 01 
O.IOOOOOOE 01 
O.IOOOOOOE’ 01 
0.2034000E 00 
0.7034000F 00 

0.0000000 
0.0000000 
o.oaooooo 
0,0000000 
•0.1033400E 00 
OrTOT^4-Ot)E 00 
0.1298770E 00 

0.1298770t^ 00 

0.9276694E-01 

0.9276694E-01 

0,3^73000E-01 

0-s-5’3-7-*?:0 OOE -01 - 

0.5375000E-01 

0.5375000E-01 

0-.rT4550t)E-01- 

0.1745500E-01 


'0'."57437&3fE=tir 

O.OOOCCOO 

0.00^0000 

o.ooooooo 

0.0000000 

-0“.1O0cjCO0-E“O1 
-0.2438S16E 00 
0.5093200E 00 
-0.9595554E--0Z 
0.3983499E 00 
-0. 2246330E-02 
0.'2914700E' 00 
-0. 136^395E-02 
0.2171200E 00 
-n,70] 1979E'-0'^ 
0.1641000F 00 
0.5580829E-03 
0.-r64400CTE“00 
-0. 1022857E-05 
0.4040000E 00 
0. 166609 8-E^02 
O.IOOOOOOE 01 
0.3970001E-02 
O'i rooooooE~oi 
0.0000000 
0.2034000E 00 
-0.1229886F-01 
0.0000000 
-0.1253360E-02 

O.OOOuOOO 
O.OOOOOOO 
-0.'10334CrcrE“0'0 
-0.2119718E-02 
-0.1298770E 00 
- 0 .-4*'4 51-0 3 0 3 

-0.9276694E-01 
-0.1495S25F-03 
-0'.5^75000E-^1 
-0.764321 3E-03 
-0.5375000E-01 
--OTlli-l-tO^TrE^t)^ 
-0. 1745500E-01 
0.4995491E-03 
O'. 48tl330'OE=01 
0.3093649E-03 


-0 .‘5'67roT8FF-0'5 
0.1072000E 04 
0.0000000 
'0.11 6-5000 E '04 
0.0000000 

OT262trOOO'E 01 
0.1731974E-01 
0.7250000E 01 
0.rS21215'E-04 
0.3500000E C2 
0.8744358E-05 
0.-6500000F 02 
0.3927955E-05 
0.1150000E 03 
-0.'5640'276'E-06* 
0.1650000E 03 
-0.9906965E-05 
0-;23'0‘DOOOF '03 
0.1876347E-08 
0.7055000E 03 
0.-3891848E-05 
0.7650000E 03 
-0.6600013E-05 
"0T'86'50Xl0O*E 03 " 
0,0000000 
0.1029000E 04 
cr.-6324“489'E-04 
0.1092000E 04^ 
0.2484331E-05 

0.2620000E 01 
O.OOOOOOO 
■OT72 5'O'O'D0E 01 " 
0.6191885E-05 
0.3500000E 02 
-0TT2'9'3 4"5 2F-- CT6‘ 
0.6500000E 02 
0.6942437E-05 
"OT11-5'000-OE“0? 
0.1539321E-04 
0.1650000F 03 
•0'ir2441‘40FE=12~ 
0.6650000E 03 
-0.2120494E-05 
■■0176513'aO'OF' 03" 
-0.4855502E-06 


0.1165000E 04 
0.r293300E 04 

ff.T2 5Cn01D'OE 01 
0.3500000E 02 
0.6500000E 02 
o;ii5oocroE 03 
0.1650000E 03 
0.2300000E 03 
0.705500’OE' 03 
0.7650000E 03 
0.8650000E 03 
'0TT079O0T)E '04 
0.1092000E 04 
0.1293300E 04 

0.7250000E 01 
0.35'00'0'OOE 02 
0.6500000E 02 
0.1150000E 03 
0,1,65000'0E 03’ 
0.6650000E 03 
0.7650000E 03 
O.'OSTOOO'OE 03 
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TABLE I (Continued) 


o ^ 

4 

0,4803300F-01 

0.0000000 

0.9650000E 03 

0.1072000E 

04 

-0,40nooaO=' 

02 

0,'48'03T300F-Or 

7)T5743T63r-lT3 

-O.SIS'OOIS'EE-OS 



10 1 

4 

0.0000000 

O.OOOOOOO 

0.1072000E 04 

0.1165000E 

04 

-n 

07 

'^,0090000 

9, Of^nonno 

0.0000900 



•> 1 ^ 

4 

0,0000009 

9,0900000 

0.1T65000E 04 

0V129330OE 

04 

-A,•5i.'^oonop■ 

07 

O.OOOOOon 

0,0000000 

0.0000000 



1 3 

4 

0.0000000 

■0,-0000000 

0.2-6-20000E -01 

0-, 1-6*5 OOOOE 

'03 

'^.9f>OOOOOP 

^2 

0.0000000 

o.ooooooo 

O.OOOOOOO 



7 

4 

0,0000000 

0,0000000 

0.1650000E 03 

0.1293300E 

04 

0.9000000F 
1 0 

^2 

0.0000000 

OvOOOuOOO 

0.0000000 



1 3 

4 

0.2144500E 01 

O.IOOOOOOE 01 

0.2620000E 01 

0.7250000E 

01 

0 .ne3000E 

02 

0v2144500F- 01 

-0-.-24338J6-F“ 00 

-0»lT31r9'74E-"01 



7 ^ 

4 

O.IOOOOOOE 01 

0.8095290E 00 

0.7250000E 01 

0.3500000E 

02 

,137^900'=’ 

0? 

O.lOOOOnnc 01 

-0.9595554E-0? 

0.1321216E-04 



•a -a 

4 

0.5093200E 00 

- 0.3983499E 00 

0.3500000E 02 

0.6500000E 

02 

0 ,3900^^001^ 

0? 

9,3993200E 00 

-9. 2246330F-02 

0.8744358E-05 



4 7 

4 

0,39E74oof 00 

0,2'’14700E 00 

0.6500000E 02 

0.1150000E 

03 

0, '57‘^ 0000c- 

02 

0.-3983-49OE 00 

-0,1363-^95F^-02- 

0.39-27935E-05 



S 3 

4 

0.2914700E 00 

0.2171200E 00 

0.1150000E 03 

0, 1650000E 

03 

0.6930000F 

C2 

0.2914700E 00 

-0.7011979E-07 

-0.5640276E-05 



6 3 

4 

0.-2171200E 00 

0,l-22t'u00E 00 

O.I 6 - 5 OOUUE 03 

O'. 2300000E 

03 

0,82?5000E 

02 

0.2171200E 00 

0. 1205775E-02 

-0.1987146E-04 



7 3 

4 

0.1220000E 00 

0.7999998E-01 

0.2300000E 03 

0.265C000E 

03 

0.9600900E 

0-2 

■0,3 220000E-^ 00 

0.-298 77 58t-CT2- 

■ - 0 .-6*83 38 2 4-E-O 4 



0 3 

4 

0.7°999QeE-0] 

0.4500000E-01 

0.2650000E 03 

0.3000000E 

03 

0*]0l0C00F 

03 

0.799Q998E-01 

0.3938776E-02 

-0.8454810E-04 



9 9 

4 

o.^^onoo^E-oi 

OiOOOOOno 

O.’SOOOOOOE 03' 

0*.3650000E 

03 

''.1050000'" 

03 

0,43000noF-01 

0.7455617E-03 

-0.1119708E-04 



1 0 3 

4 

0,0009000 

0,0000000 

0.3650000E 03 

0. 1293300E 

04 

O.^O^lOCOOc 

03 

ovoe-ooooc 

'0.0000000 

OTOOOOOOtr 



1 1 
1 

4 

0,0000090 

0.0000000 

0.2620000E 01 

0.7250000E 

01 

O.ODOOOOO 


o.ooooooo 

0.0900000 

0.0000000 



7 3 

4 

3,0000000 

-0.1033400E 00 

0.7250000E 01 

0.3500000E 

02 

n,0000000 


0.0000000 

-0. 2119718E-02 

0.6191885E-05 



- "a - -c< 

- 4^' 

■-Chrt033-400€ 00' 

--0.-1-298-770E— 00 

Or3-5croooo-E 02* 

ai6300000E' 

02 

-0,1 500000E 

r»i 

-0.3033409F 00 

-9.448in35E-03 

0.1293452E-06 



i. 3 

4 

-0.1298770E 00 

-0.9276694E-01 

0.6500000E 02 

0, 1150000E 

03 

-O.FnooOOO*^ 

01 

^O-il 298770^ 00- 

■7-0. 149‘='825'E=-'03 

0.6'9424-37E-tl5' 



c 3 

4 

-0.9276694E-01 

-0.3?75000E~01 

0.1150000E 03 

0. 1650000E 

03 

-o.n oooooc 

0? 

-0,92766Q4'=--01 

-0, ■'6*^21 3E-03 

0. 1539721E-04 



6 “'3 

*4~ 

-O-S-5-37-5-O&0E-01 

'.•6hi'4999^ 9 9 F'“0‘2 

0?-l'650000e-03-- 

0v23 00000 E- 

03 

-0,1 562900C 

02 

-0.5375000E-01 

-0. 1875739E-02 

0.1795629E-04 



7 3 

4 

-0.4999999E-02 

0.9999995E-01 

0.2300000E 03 

0.3650000E 

03 

-0,2Q00000F 

02 

-0-.499999^E“02 

-0-.-4 85596-7 E«^0-3“»0. 4-T7-5 Tt9 E'- 0 6' 
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TABLE I (Concluded) 


8 ■ 3 " 4' 0.9-99g995E-0l 0''."1479999E“Oa 'T}r3’6^0000E"T)3' “O.TfS^OOOO'E 03 

-0.1300000F 02 0.9999996E-01 -0. 7799980E-03 0 . 6799986E-05 

934 0.1479999E 00 0.1059999E 00 0.4650000E 03 0.5650000E 03 

~ 0 , 4000 '' 00 '^ r>l” 0.r4"T999'5E 00 ~ 0 . 180 C 01 5E=^03 - 0 ♦ ZeTOOlOE-tTS 

■[0 3 4 0.1059999F 00 0. 000^000 0.5650000E 03 0.7650000E 03 

n .inoPCOOF 07 0.1059999E 00 0 . 6500068E-04 -0 . 1 100001 E-05 

n - 3 - 4 - O.-OHO-OOOO “O-.OOO'JCOfO “ 0 r7650000E' 03 ’0.rZ93300E 04' 

n. TKQOOOO'^ no n.oOOOOOO ^.0000000 0.0000000 

2 

1 3 4 0.0000000 O.OOOtOOO •“0V26'2'0000E 01 'C.r6'50000F 03 

O.9O0nOCOE n 2 O.OOOOOOQ O.OOOUOOO 0.0000000 

234 0.0000000 0.0000000 0.1650000E 03 0.1293300E 04 

O.QDOOOOO'^- C? ' n, 0000000 O.-ODOOCOO '0.0000000 

6 

134 0.0000000 0.0000000 0.2620000E 01 0.3500000E OF 

o. nonneoo n.ooo^'Onn n.tiOOOOOO "O.OOOOOCO 

334 0.0000000 0.1603999E 00 0.3500000E 02 0.6500000E 02 

0, 3003000 0.0000003 0.598C0O0E-02 -0.7340743E-04 

-3 4 -'0.10O3OOOF- 00 — 0','3a-750aor"00~ -O.'SEOOOOOF 02 O.ieEOrOOOE 03 

0.3400000“^ 01 0.160K9OQC- 00 0 . «930074E-03 0 . 1099 B45 E-06 

i 3 4 0.342'iOOOF 00 0. 34280002 00 0.1650000E 03 0.3650000E 03 

0.7830000P 02 O. 34 - 75 OOOE 00 -0vn44408E=O8 ■-0.58r4'695E-ri 

334 0.3425000E 00 0.0000000 0.3650000E 03 0.4350000E 03 

0 , <^700000'^ 02 G. 3425 OOOE 00 -0.3663260E-02 0 . 1 1 58888 E-04 

^ - 4 - O.OO&OOOO -0.00000r)0 ' -0."4-35ff00aE -03- OTl-2933'OOE 04 

0.1070H00E n.3 0.0000000 , O.OOOUOno 0.0000000 

10 

1 -3- 4 0.2r4'4500E 01 0. 96569O0E -QO 0.-2620000E 01 0.7250000E 01 

O.1183000<^ 02 0.2144500E 01 -0.2014851E 00 0. 1068164E-01 

7 3 i O.9656900E 00 0.50°5299F 00 0.7250000E 01 0.350C000E 02 

0.1030COOE n-> 0.'9fr8'S900E' 00-^0."3T2'72E3E'-02~=t).lT99?35^E-03 

384 0.5095299E 00 0.3983499E 00 0.3500000E 02 0.6500000E 02 

n,4n7«;000F 02 0.509529Q'=’ 00 0.7752996E-02 -0.2134667E-03 

4 3 4' 0,39&3'499E'00 ' 0.-2914T00E“00'" O.'&SD'OaOOE 02 O'. lT5‘000"0r 03 

0.5675000F 02 0.3983499E 00 0 . 5236595E-02 -0.8407197E-04 

534 0.2914700E 00 0.2171200E 00 0.1150000E 03 0.1650000E 03 

O.797 5OO-0«^O2' O', 2'^rKTOOE~'Ot) ■0T823 3TS'9E^2 ~=tmT9T639E-03‘ 

634 0.2171200F OO 0.3939100E OO 0.1650000E 03 0.1950000E 03 

0.9945000E 92 0.2171200E 00 0. 756 1672E-02 -0 . 1025595 E-03 

7 8- - i 0, "3-g 3 - 9 - 1 OOE' 00 0','9489F50E~00 '~07T93'0000E 03 0,1993000'E '03 

O.UOOOOQF 03 0.3989100E 00 0.2066403E 00 -0.2147668E-01 

834 0.9489650F 00 0.9489650E 00 0.1995000E 03 0.23450C0E 03 

0-i 1 1-40000F- 03 0-^, 94 8-963 0-E“00--=Ot5E3^4“3E=^3 -0t9®'7T3 2 4'E='03' 

934 0.9489650E OO O.OOOOOOO 0.2345000E 03 0.2790000E 03 

0. 14700008 03 0.9489650F 00 -0 . 1841072E-01 0, 1160779E-03 
10 "'8-' - 4- OiOOa&OOO — O.-tlOOOOOO -0T2T9tr00t>E~03' '0;'T293300E 04' 

0.1630000F 03 O.OOOOOOO O.OOOsjOOO O.OOOOOOO 
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TABLE II 


Labels for logical unit numbers used for input and output of data 


from disk 

or tape 

files : 


Proqram 

1 

(FVTM) 



Label 



Unit No, 

Purpose 

INPUT 



5 

input from punched cards 

OUTPUT 



6 

output on printer 

GEOM 



9 

body geometry data 

P ID AT A 



10 

storage output 

PlIN 



11 

restart input 

Proqram 

2 

(FVSM) 



INPUT 



5 

input from punched cards 

OUTPUT 



6 

output on printer 

GEOM 



9 

body geometry data 

P ID AT A 



10 

input of program 1 storage file 

P2RSIN 



11 

restart input 

P2IFIN 



12 

unfinished interface storage 
file input 

P2RS 



13 

restart output 

P2DATA 



14 

interface storage file output 
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TABLE II (Concluded) 


Procrram 

3 (-PDSM) 



Label 


Unit No. 

Purpose 

INPUT 


5 

input from punched cards 

OUTPUT 


6 

output on printer 

DATA20 


20 

input or output on dish 
(binary) 

DATA21 


21 

input or output on disk 
(binary) 

CONTIN 


22 

input from continuous storage 
file on disk 

CONTOU 


23 

output of continuous storage 
file on disk 

PUNCH 


30 

punched card output 

BBDATA 


31 

input from program 5 
(blunt body data) 

Procrram 

5 (CTI) 



INPUT 


5 

input from punched cards 

OUTPUT 


6 

output on printer 

GEOM 


9 

body geometry data 

PIDATA 


10 

input of program 1 storage 
file 

CHAOS 


11 

output for program 4 (CHAOS) 

BBDATA 


12 

output for program 3 (FDSM) 

P2DATA 


14 

input of program 2 interface 


storage file 
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TABLE III 


JOB CONTROL CARDS TO RUN PROGRAM 1 
FROM AN INITIAL START 


Card 

1 

2 

3 

4 

5 

6 

7 

8 


JOBNAME, T , P_, 

ACCOUNT, ID NAME, ACCOUNT #. 

ATTACH, LGO, PROGRAMl, ID=AAAAAA, PW=BBBBBBBB. 
ATTACH, GEOM, SHUTGEOM, ID=AAAAAA, PW-BBBBBBBB. 
REQUEST, PIDATA, *PF. 

LGO 

EXIT, U. 

CATALOG, PIDATA, DATAl, 1D=AAAAAA, MR=1, RP=120. 



6 



INPUT FOR PROGRAMl (initial start) 
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TABLE IV 



JOB CONTROL CARDS TO RESTART PROGRAM 1 AND RUN 
ALL SHUTTLE FLOW FIELD PROGRAMS 

Card 



1 

JOBNAME, T , P_. 


2 

ACCOUNT, ID NAME, ACCOUNT #. 


3 

ATTACH, LGO, PROGRAMl, ID=AAAAAA, 

PW=BBBBBBBB. 

4 

ATTACH, GEOM, SHUTGEOM, ID=Ai\?JiAA, 

, PW=BBBBBBBB. 

5 

ATTACH, PlIN, DATAl, ID=AAAA7^. 


6 

REQUEST, PlDATA, *PF. 


7 

LGO. 


8 

EXIT, U. 


9 

CATALOG, PlDATA, DATAl, ID=AAAAAA 

, MR=1, RP=120 

10 

RETURN, LGO. 


11 

RETURN, PUN 


12 

ATTACH, LGO, PROGRAM2, ID=AAAAAA, 

PW=BBBBBBBB. 

13 

REQUEST, P2RS, *PF. 


14 

REQUEST, P2DATA, *PF. 


15 

LGO. 


16 

EXIT, U. 


17 

CATALOG, P2RS, DATA2, ID=AAAAAA, im=l, RP=120. 

18 

CATALOG, P2DATA, DATA3, ID=AAAAAA 

, MR=1, RP=120 

19 

RETURN, LGO. 


20 

ATTACH, LGO, PROGRAMS, ID=AAAAAA, 

PW=BBBBBBBB. 

21 

REQUEST, CHAOS, *PF. 


22 

REQUEST, BBDATA, *PF 


23 

LGO, 
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TABLE IV (Concluded) 


24 

EXIT, U. 





25 

CATALOG, CHAOS, DATA4, 

II 

Q 

H 

AAAAAA, 

MR=1, 

RP=120 

26 

CATALOG, BBDATA, DATA5, 

ID 

=AAAAAA, 

MR=1, 

RP=120. 

27 

RETURN, LGO. 





28 

ATTACH, LGO, PROGRAMS, 

ID= 

AAAAAA, 

PW=BBBBBBBB. 

29 

REWIND, BBDATA. 





30 

LGO. 





31 

EXIT, U. 





32 

REQUEST, CONTOU, *PF. 





33 

CATALOG, CONTOU, DATA6, 

ID 

=AAAAAA, 

ME=1, 

RP=120. 



INPUT FOR PROGRAMl (restart) 



INPUT FOR PROGRAM2 



INPUT FOR PROGRAMS 



INPUT FOR PROGRAMS 
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TABLE V 


JOB CONTROL CARDS TO RESTART PROGRAM 2 AND 
RUN REMAINING SHUTTLE PLOW FIELD PROGRAMS 


Card 


1 

JOBHAME, T , P_. 



2 

ACCOUNT, ID NAME, ACCOUNT #. 


3 

ATTACH, LGO, PROGRAM2, 

ID=AAAAAA, 

PW=BBBBBBBB. 

4 

ATTACH, GEOM, SHUTGEOM, 

, ID=AAAAAA 

, PW=BBBBBBBB. 

5 

ATTACH, P2RSIN, DATA2, 

ID=AAAAAA. 


6 

ATTACH, P21FIN, DATA3 , 

ID=AAAAAA. 


7 

REQUEST, P2DATA, *PF. 



8 

LGO. 



9 

EXIT, U, 



10 

CATALOG, P2DATA, DATA3. 

, ID=AAAAAA, 

, MR=1, RP=120 

11 

RETURN, LGO. 



12 

ATTACH, LGO, PROGRAMS, 

ID=AAAAAA, 

PW=BBBBBBBB, 

13 

ATTACH, PIDATA, DATAl, 

ID=AAAAAA. 


14 

REQUEST, CHAOS, *PF. 



15 

REQUEST, BBDATA, *PF 



16 

LGO. 



17 

EXIT, U 



18 

CATALOG, CHAOS, DATA4, 

ID=AAAAAA, 

MR=1, RP=120. 

19 

CATALOG, BBDATA, DATA5, 

ID=AAAAAA, 

MR=1, RP=120. 

20 

EXIT, U. 



21 

RETURN, LGO. 



22 

ATTACH, LGO, PROGRAMS, 

ID=AAAAAA, 

PW=BBBBBBBB. 

23 

REWIND, BBDATA. 



24 

LGO. 
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TABLE V (Concluded) 


Card 




25 

EXIT, U. 



26 

REQUEST, 

CONTOU, 

*PF. 

27 

CATALOG, 

CONTOU, 

DATA6, ID=AAAAAA, MR=1, RP=120. 



INPUT FOR PR0GRAM2 (restart) 



INPUT FOR PROGRAMS 



INPUT FOR PROGRAMS 
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TABLE VI 


JOB CONTROL CARDS TO RESTART 
PROGRAM 3 


Card 


1 

JOBNAME, T , P_. 



2 

ACCOUNT, ID NAME, ACCOUNT #. 


3 

ATTACH, LGO, PROGRAMS, 

ID=AAAAAA, 

PW=BBBBBBBB. 

4 

ATTACH, CONTIN, DATA6, 

ID=AAAAAA. 


5 

LGO. 



6 

EXIT, U. 



7 

REQUEST, CONTOU, *PF 



8 

CATALOG, CONTOU, DATA6, 

ID=AAAAAA, 

MR=1, RP=120. 


INPUT FOR PROGRAMS 
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TABLE VII 

JOB CONTROL CARDS TO RUN PROGRAM 4 (CHAOS) FROM A 
DATA SET PRODUCED BY PROGRAM 5 (CTI) 


Card 


1 

JOBNAME, 

T , P_. 

2 

ACCOUNT, 

ID NAME, ACCOUNT #. 

3 

ATTACH, 

LGO, CHAOS, ID=AAAAAA, PW=BBBBBBBB 

4 

ATTACH, 

TAPES, DATA4, ID=AAAAAA. 

5 

LGO. 




INPUT FOR CHAOS 



43 




Figure 1. - Mesh geometry for program 1 (fvtm) 
determined by series of nested cones. 
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Figure 2 - Partitioning of the shock layer 
into finite volumes. 






Fxgure 3,- Illustration of position, orientation, and 
vertex angle of conical data surfaces for 
successive steps (program 2, FVSM) . 
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Card ETo 1 LOGICAL VARIABLES AND INTEGERS 


Variable 

Card Column No 
Data 


Poxinat UL1»3X»I5,X1) 


B 


[ 


nw. 


HIMHEI 

1 


II 


■Bl 

1 



Card No. 2 INTEGERS-RIGHT ADJUSTED - 5 column fields 


Variable 
Card Colxunn No 

Data 


Format ^12IB) 


NPRT 

NEND 

IE 

JE 




5 

10 

IS 

20 

25 


IS. 









Card No 3 REAL NUMBERS - 10 column fields^ decim^al point required 






Variable I 

Card Column No 
Data 1 




GAMMA 

SDO ' 

C2 1 

g? — 

10 

20 

20 

40 

50 

60' 










REAL NUMBERS - 10 column fields, decimal point required 


Variable 
Card Column No 



7.FOCNL 

THWL 

THLL 

RBX 

10 

20 

30 

40 

50 







(a) Program 1 (PVUBB) 
Figure 6 • Input data form 






O) 

o 


card Wq — 1 imtEGERS - I colisnn fields 
Format (5I1> 


Vara, able 


Card Column No 
Data 


. NRSTRT 

NPfLOW 

— NPAREA 
rp-NPXY2 
, J I ^ NPPIUX 


EEUas 


”° — I REAI- number - 10 colunn field, 

Variable 
Card Column No 
Data 


Format {Flo 0) 


dthetp 
10 


decimal point required 


(b) Program 2 (PVSbb) 
Figure 6 - Continued 




Card 1 LOGICAL VARIABLES, INTEGERS (RIGHT ADJUSTED) AND REAL NUMBERS 


Variable 

Card Column Ho 

Data 



1 MODIH 

1 l/rlN 

1 HODOUT j 


■eai 

HODCPO 

■JW0=H.l 

CRASHZ 






iiiH 




1 1 










card NO 


INTEGBR-RIGJiT ADJUSTED 


Format (15) 


Variable 

card Column No 

Data 


NUMDUM 


5 column field 


Card Nos 3 to 


Variable 

Card Column No 

Data 


(NUMDUM + 2) REAL NUMBERS - 10 column fields, decimal point required 
Format (7510 01 


ZRMSH 

BETTA 

PHIZRO 

'AR 

BR 

ZPACT 


10 

20 

30 

40 

50 

60 










Card No (NUMDUM 4- 3 } INTEGER-RIGHT ADJUSTED - 5 column field 


Variable 

Card Column No 

Data 


Format 

NMPTS 

5 


(15) 


(c) Program 3 
Figure 6 - Continued 







card Nos 


<MUMDUM4-4) to (NUMDUM-f tJHPTS 4- 4) REAL NUMBERS, 
^ V required 

Formab (2ElO 0) 


Variable 

Card Column Mo 

Data 


ZDMPM 

DCM 

10 

20 




10 column fields, decimal point 


Card No (MUMDUM -h NMPIS -f 5) INTEGER-RIGHT ADJUSTED - 5 column field 
Format (15) 


Variable 

Card Column No 

Data 


NRPTS 
5 


Card Nos (NUMDUM + NHPT5 + 6) to (NUHDUM+ NMPTS 4- NRPTS + 6) REAL NUMBERS 
10 column fields, decimal point required 


Variable 

Card Column No 

Data 


Format (2E10 0) 


ZDMPR 

DCR 

10 

20 




(c) Program 3 - Concluded 


Figure 6 - Continued 




Cana tto 1 INTEGERS-RIGHT ApjuSTED - 5 column fields 


Variable 
Card Column tfo 

Da Ca 


Format (1015) 


JLIMIT 

JL2 

NUMK 

lOPT 

NPRNT 

S 

10 

15 

20 

25 







Card Ho 2 INTEGER-RIGHT ADJUSTED - 5 column field 


Variable 

Card Column No 

Data 


Format (15) 

12 

5 


Card Nos — 3 bo 7 as needed f or up to 40 values REAL NUMBERS - 10 column fialds, decimal point required 


Format (8Fl0 0) 


Variable 

Card Column No 

Data 

ZM 

rRC 

1 

10 

20 

30 

4(2 

50 

AO 

70 

80 










































(d) Program 5 (GTl) 


Figure 6 - Continued 






Card No 1 REAL NUMBERS “ 15 column fields, decimal point required 
Format {2B15 7) 

Variable 
Card Column No 
Data 


RNOSE 

CONAHG 

15 

30 




Card No 2 


INTEGER-RIGHT ADJUSTED 
Format (IS) 


Variable 
Card Column No 


NSC 


5 column field 


Card No 3 INTEGERS AND REAL NUMBERS 


Variable 

Card Coliunn No 

Data 


Format (3I5,4 e15 7) 


NS 

K 

NCF 

YPl 

VP2 

ZL(NS.H. 11 

. N. 21 

5 

10 

15 

22. 

45 

GO 

75 









Card Nos 4 to 6 REAL NUMBERS - 15 column fields, decimal point required 
Format (5E15 7) 


Variable 

CF 

I 

Card Column No 

IS. 

.30 

45 



Data 













(e) Body geometry input data 


Figure 6 - Concluded 






Figure 7.- Body and shoclc geometry rn pitch plane ‘Showing 
definition, of computational mesh input parameters. 




Figure 8. — Estimated axial iQcation of last 
mesh 'cone vertex for program 1 
(see figure 7) . 





Figure 9.- Control points and geometric parameters for the body geometry 









FT 01 


1 

200 16 12 

18 10 -1 




10 0 

797.7905 1 

0269E-06 1.4 

.26 

.5 

o 

• 

o 

30.0 

7.0 

40. 90. 

1,0 




(a) 

Input data for 

program 1, 

(FVUBB) 


01111 






5. 







(b) 

Input data for 

program 2. 

(fvs’bb) 



11 

22 17 

43 1 






40 








.025 

.05 

. 075 

.10 

. 125 

.15 

.175 

.20 

. 225 

.25 

.275 

. 30 

. 325 

.35 

.375 

.40 

.425 

.45 

.475 

. 50 

. 525 

.55 

.575 

.60 

.625 

.65 

.675 

.70 

.725 

.75 

.775 

.80 

.825 

.85 

. 875 

.90 

.925 

.95 

.975 

1.0 



{c) Input 

data for 

program 5 

(CTI) 




FALSE 

5 

1 

3 

0 0 

50 

200. 

5 

10. 


5.0E. 2 

O.OE.O O.OE.O 

0. OE.O 

0. OE.O 

5.3E.2 

0 OE.O O.OE.O 

0. OE. 0 

1. IE. 0 

5.7E.2 

3.0E.01. 30899E.0 

O.OE.O 

I.IE.O 

6.0E.2 

3.0E.01.30889E.0 

0 OE.O 

1. lE.O 

3.0E.3 

5.0E.01.30899E.0 

0 OE.O 

1. IE. 0 

2. 5E. 2 

0. OE. 0 



3.0E.3 

2 

4. OE-1 



4. 5E. 2 

O.OE.O 



3.0E 3 

1. OE-1 




(d) Input data 

for program 

3 (FDSC) 


Figure 11.- Input data for niomerical example. 


60 



014 

00008017 
003002 
020010 
000000 
000000 
000000000 
107.3 90 

40 2 2 333 

1.4 1.4 1716. 

1. 3639E-011.9701E-07 21700. 


89.056 

0.0 

0 

180. 

0. 0 

0.0 

179.8 

0 0 

0.0 

179.6 

0.0 

0.0 

179.4 

0.0 

0.0 

179. 

0.0 

0 0 

178.5 

0. 0 

0.0 

178 

0. 0 

0.0 

177.5 

0. 0 

0.0 

177. 

0.0 

0.0 

176. 

0.0 

0.0 

175. 

0. 0 

0.0 

174. 

0. 0 

0.0 

172. 

0 0 

0.0 

170 

0. 0 

0. 0 

168. 

0. 0 

0.0 

165. 

0. 0 

0.0 

160. 

0. 0 

0.0 

150. 

0.0 

0.0 

130. 

0.0 

0.0 

0.0 

0.0 

0.0 

SHUTTLE 

147 


PERFECT GAS 



000 



000 



. 0726 

109 

. 145 


0.0 1.5 


.182 .218 


(e) Input data for program 


. 254 


(CHAOS) 


Figure 11,- Concluded 



3"0|HfN!)t0NAl UHlfjr nnpY tL 

MACMil .oonnnnt (PI '■pps*? PTriosi^o? unosi ,(pSft'>p{it>.ti6 CAXHAst.iinggoottno 

SOo*2 <<ognPn(>n! CPI's nnAongf.ni 'sjm, Au«si(goggooF»gi Afcc>ii.>F|0ooo(>0Ea<iO 

I HCNO. ?n Ifl Ih JH 12 Kt« 16 JSHKi 10 *.ei5i .1 
THW| s A0606pC*0? IHU I ,9(loggol<Op «BX I (lOOOOnEAOl 

VINF« 3 2?Ti>'iF«rii; AlMFi i ?1i65r»fls H!6Fi p,Tlgl2t«09 BM/i 3 , 366«Tt-l)t tlNfi S,Ta396tA06 EHNF9 |,91I223E«0« Tl«F» 2,H567E*0J 


npOCCK.I.n CON3TA9T VALI'F? 

1 ? F 

I? 13 7“ 

10. 0 .1 o*‘>r40« 

-I,OS«E»OP -t 3l7f4(l2 .l.lpaF.O? 

2 0, 6, .1 nsOE+OO 

-I.oooeaOp -i,3t7r*o? .i.SkaF.si? 

3 0, 0, .t,0«<7E*l)O 

• l,«aoE*0P •I,3|7F«H2 .1 Sf.8F.02 


I) 

5 

6 

T 

e 

9 

10 

11 

12 

13 

10 

IS 


0. 0 .l.OftSF.OO 
■I.OOOE.OP -1,X17F*02 .I.SfSE.O? 
0. 0, .1.089F.0O 
il.oSOE.O? -l,SI7F.(l? .l,Sf,6£.02 

0. 0, .I.OaOL.OO 
'i.oeoE.oa »i,3i7F*o? .i.Saof.o? 
0. 0 .I,0s9f.0l> 
il.089t.0P -1,317F.0? .1,5 a8E«OS 
0. 0 .I,Oa 9C«OI7 
.1.089E.0? .l,ll7F.n? .l,5*nF.02 

0. 0. .1 OpSE.IlO 
il.OB9E.OP .l,317F+0J .I.SpBE.OJ 
0. 0, ,I,Oa9E.OO 
1.089E.O? .1,3I7F.0S .1 SaBE.02 
0. 0, .1 OB9E.OO 
I.O09E.OP >l,tl7F.02 .l.SpBE.OF 
0. 0. .1.O89E.00 
1.009E.05 -IfllFEtOE ,1 ,5 a8E.02 
0. 0 .l,nA9F.oo 
1.089E.OP -l,Ft7F.O? .1 S(>nF*li? 
0. 0( .1 OpOE.OO 
1.069E.0? -1,A|7F*0? .l,S(,eF*02 
0.- 0, . .1 OR9E.OO 
I.089E.0P i.l,3l7E.oa .1 SRBE.OP 


. 10 

17 

18 


0. — 0, .1.0a 9E«Q0 
1.089E.0? .1.3I7F.02 .1,5*BF.02 
0. 0, .1 089C.O0 
l,089E.0p -1,317F.02 .I,Sa 8E.O? 
0. 0. .1.089F.00 
l.OBoe.OP -1,1I7F.02 .I.SaSI.OS 


BF0C(K,2,I) C0NS71NT VALKfS 

1 2 3 

12 . 13 i« 

1 T.miE.Oj 7,M9F.0| 7,3n3E.01 

2.597E.0? 2,903F.|)2 5,‘lp6E*02 

J I.lUE.Ol 7,MaF.ni 7 30SE.01 
2.S97E70? 2.993F.02 3 0 |,ae.o? 

} 3>|10EAOi 7 MOF.OI 7,3o3E.0l 

2.S97E.02 2,993F*02 3 Oi6e.02 

A 7.1UE.0J 7,tlOF.<l| 7(3o3E.01 
2.S97E.02 2 995F.I12 3,aihE*02 

* 7,)UE.0| 7,(l«F.0l 7 3(i3t.0t 

2,S97t.0j 2,»93F«02 l.OlOE+O? 

0 7,i19E.0i 7,||ilF*ni 7,3(i3F*01 

2.S97E.02 2.993F.02 3 «p6t.02 

7 7.IIOE.OI 7 tlUF.Ol 7,3o3E.O| 

2.597E.0f 2.993F.02 3 OpOE.oa 

8 7,414E*0i 7,(I<IF.Oi 7 SgSE.Ol 

2.S97E.02 2.993F.02 I.OptE.OP 


a 3S9E(0(i .9,797b>o9 -i,7U2t.nt -PiFaiE.gl -3.9I9F.01 i.S,339t.0l >9,9072*01 

1 Boot.op 1392.(12 

il,3S9E(00 .9(797t.o0 >l,7U2f.oi >2,721E.0I >3,9192.(11 .S,33OE«01 .0.96TE.01 

l,89(iE,02 >2,1391*02 

9.3591*00 >9,7972*00 >t,792E.01 >2,7212.01 •5,919E*01 >S,339E*01 >6.9li7E*0l 

l,«90t,02 -?,t3“E*o2 

9 354E,60 .9,7972*00 *1.7922*01 >2,7212*01 >3,919E*0I >S. 3392*01 .fc.967E+0l 

1,890E.02 >2.1392*02 

0. 359E,00 .9, 7972*00 *1,792E*0I >2,72|E*0I >3,9l9E*01 >S,339E*01 >6,9(,7f*0t 

1,80oE,02 .2,1392*02 

*’.7972*00 >1,792E*01 >2,7212*01 >3,919E*01 >5,3392*01 >9,967F*0I 

1. eooEloa >2,139E*o2 , ' 


II 

>8,8172*01 

>8,B17E*0I 

>8,0172*01 

•8,0172*01 

.B,8|7E*0I 

>8,0172*01 


iO,3S9E*00 >9,797E*nO >1.7922*01 -2,7212*01 >3,9|9F*0I -5,5392*01 >6,987E*01 *6, 0172*01 
■1,890F,02 .2,1392*02 

.9.3S9E*00 -9,7972*00 >1,7922*01 -2,7212*01 >3,9192*01 >5.5392*01 >6,9672*01 -8,8172*01 
it. 8902*02 -2,1392*02 » < ‘ 

■9,3S9F*00 >9,7972*00 -1,7022*01 >2,7212*01 -S.OlOE.Ot >5.3302*01 >6.967F*0l •0,0171*01 
il,890E,02 >2,1392*0? . 

■9,3S9F,00 .9.7972*00 >1,7922+01 >2,7212*01 •!, 9192*01 -5.1392*01 >6,9672*01 -8.8172*01 
■1, 8902102 >?, 1392*02 ' 

.9.3592,00 >9.797E*00 >1,792E*01 -2,7212*01 -3,9192*01 >5,3392*01 -6,»672*0l .6,8172*01 
•1.8902102 -2.1392*02 

•9,3592*00 >9,7972*00 >1,7922*01 >2,7212*01 •1.9t9£»ol -5,339E*0I >6,9672*01 >8,8172*01 
•1,890E1 o 2 >2.1592+02 i 

9.3592.00 >9,7972*00 >1,7922*01 >?,72IE*01 >3,9192*01 >5,1392*01 >6,9672*01 >8,8172*01 
1,8902*02 .2.1392*02 

9.3592.00 .9,7972*00 >1,7922*01 >2,7212*01 >l,919t*0l >5,3392*01 >6,9672*01 >8,6172*01 

•1,8902102 .2. 1392 + 02 i • ‘ 

•9,359E,00 >9,7972*00 >1,7922*01 >2,7212*01 >3,9192*01 >5,3392*01 >6,967£*01 >8,8172*01 

1.B9o2;o2 >2,1392*02 

9,3S9£j00 .9,7972*00 .1,7922*01 >2,7212*01 -1, 9192*01 >5,3392*01 -6,9671*01 >8, *172*01 
1 BOOtJoE .2,1392*02 ' 


9,3592,00 .9,7972*00 -1,7922*01 -2,7212*01 >3,919E*0t >5,3392*01 .6,9672*01 >8,ei7E*0l 
1,B90E^02 >2, 1392*02 

9,3592*00 .9, 7972*00 >1,7922*01 >2,7212*01 >3,9|9E*0I -5,3392*01 .6,9672*01 -8,8172*01 
1,8902,02 .2,1392*02 


4 

5 

6 

7 

IS 

18 




8.8112*01 

l.OUE+02 

1.183E*0? 

3.898E«02 

9,9072*02 




8,8112*01 


U1S3E«0Z 

3 «’>8e4.02 

9,9072*02 



7,a6«E40i 

8,8UE*ol 

l.JH3E4r02 


3,B9aE*02 

0.907£*o2 



7 ^ 8as£^o ] 

8 . 8112*01 

1.013E*02 

(■IS3E*0£ 

3,B08F*02 

0,91)72*02 



7*8<»8£;0i 

8 , 8112*01 

U9t3Et02 

tiiesE-ftOZ 

5,e^«E*02 

9.9o7£*02 




8 , 8112*01 



3,a9eEio2 

9.907F+02 




8.8112*01 

l.OlJE+02 

t 1B3E«03 

3.69dt«02 

9,9072*02 



7,8h6E^O| 

8.8112*01 

l.U1i£4i02 

Uf R3£*0? 

3,69flt*e? 

9,9072*02 




a 

9 

10 

M 

U390C402 

l,OJ5E*Oa 

1,910C402 

2,239E402 

1,3901402 


1.910E402 

2.239L402 

1«390E402 

1,635E402 

1,9|0E40£ 

2.239E402 

1«390E402 

I.635E402 

1,9i6E402 

2,239C40I 

1 ,39QE402 

1 .63SE402 

1.910E402 

2,259E40£ 

1«390E402 

1,63SC402 

1.916E402 

2,239E4Q2 

t,39l)E402 

1«63bE4Q2 

U9I6C402 

2,239E*02 

1.390E402 

t.O3SE402 

1«91A£402 

2.239E402 


Figure 12 - Selected sample output data from program 1 



9 

7,11^E f0i 

7,iiuF40! 

7 

JftU *01 

7 


1 ,Bt3t*02 

l,]83F*02 

1 ,300£702 

l,6iSE *0? 

l,*»lBf 40a 

2,2i«t402 


2 •»97f*0? 

j ao^F*Ai» 

\ 

*0? 

3 80>»E*it2 

■ ,«07U(i? 






10 

7,ilxie»lH 

7 « 1 1 ur*m 

7 

Sft3f *.M 

7 *f»i 


1 .iJt Jr*rt? 


1 ,3O0f *02 

t,b3^F«0P 

),RtBE402 

2,2J9fc*0a 


2,S97f*0p 

? coj^ + n"* 


*0? 

3,8‘>flE,02 

q,al)?Ho? 







n 

7,1 let fOl 

7 iiHF+ftl 

7 

SftSl *fl1 

7 ,Oj 

» «m»oi 

l,(il3E*02 

1 16SE*02 

1 ,sooe*o2 

l,63St*02 

l,9lAf*02 

2,2S9E402 


2,577P*02 

2,oo^F+r>? 

s 

«?fet+02 

i,8'’^E*02 

•o^ 







12 

7, 1 iaF40» 

7 lMF+fH 

7 

SrtSf *C1 

7 

»,i'iif*(ii 

1.0I3E*02 

1 , |H<E*0P 

],SO0E*f>2 

1 ,63Sfc*0? 

|,<*iab*02 

2,219E40^ 


2*S97t*0? 

? QO^F,fi? 

% 

4^6^402 

^ fl«8E*Ci2 







IS 

7*liaE4Di 

7 IMP+Ol 

7 

Jo'f *0J 

7,ftortF,Oi 

4,ei it«pi 

1 ,OISE*02 

1 , lB3i*02 

1,340F*02 

l,63*tfe*02 

l*Rl6B*02 

2,2i‘>fc40? 


2,*i97E402 

2,on^r*rt? 

3 

Apepf*02 

3,8R8t,02 

q,a(i7E*(i? 







lU 

7,1 |uE«0i 

7,1 WJF-flH 

7 


7,Bhflfc*0| 

4,Bllfc*(l| 

1,0|3C*Q2 

1,1BSE*02 

1,S40E402 

J,63'iE+02 

1,QtflE402 

2,239E«02 


2,«\'>7E*0? 

? 


a?6F40F 

3,HOHt+02 

o,ao?E*o? 






IS 

7,iliiE40l 

7,1 1uF4ftl 

7 

SrtSF *01 

7 ,ABbE^0t 

ii,niiF*(M 

l,0lSt*l)2 

i«iesE*o? 

t,3O0€*o2 

1 ,b3se*02 

|,^UE402 

2,2S<»E402 


2«S‘>7E*0? 

2 oujF+ft? 

S 

4?6F«02 

1 BRhfc,02 

J.4(i7E*(>2 







lt> 

7,MtfE40l 

7 iiuF+ni 

7 

SoSf*0J 

7,d6eE«0$ 

8,niit*(ii 

l,Ul3t*i>2 

1*1B5E*02 

l,590fe402 

|.63St*02 

J*4t6E+03 

2,23<>E402 

17 

2,S97E40? 

? OO-^F + flj 

3 

4?6e*o; 

i,8«eF*02 

q,uo7t»o2 







7,it^E40l 

7^n«F+fli 

7 

Sfvif *01 

7,Hf.aE40l 

e,Biu*o1 

>,0l3et02 

I tt8S£*0P 

US90E4O2 

|,t>iSL402 

1,R18E402 

2.239E40Z 


2,«>97E40? 

2,0'»S*^*02 



,02 







18 

7,1UiE40t 

7^11«F+oi 

7 

SfiSJ *01 

7,Bb8E*ni 

e,fliu*o) 

|,013E*02 

i,ias£»o? 

1,3<70E*02 

1 ,a3S£«0? 

l,9jet402 

2,239t*02 


a,*>97F*0? 

2 9OJF+02 

S 

A?hr*Q2 

*02 









RBDDY(K»T) CnNSTiNT V41 



1 

2 

3 

a 

5 

6 

7 

8 

9 

10 

U 


1 ? 

IS 


15 

19 







1 

7,]UjE«I>| 

7,1 12F401 

7 3l9E4fll 

7,9bb£*0l 

4,<l63t40l 

1,0906+02 

1 ,2546+0? 

1,4876*02 

1 T57E+02 

2.0646+02 

2,4066+02 


B.7S3t*07 

S,lO 2 F *02 

3 6aU 40? 

4,136E*02 

a, 6786+02 






2 

T.tUEtOi 

7,M2F4ni 

7,3i9t-*01 

7,9b0t*0| 

9,OS7E+ot 

l,o5fl£402 

1,2516+02 

1,462E*02 

1 7486+02 

2.0486+02 

2.3796*02 


B,7JOF*0? 

S lpaF402 

3 S<3F*0? 

a.02?E*02 

a,S276+p? 







3 

7,il<iE*l>i 

7,U?F*01 

7,5«ef+0l 

/,9SB£*oi 

9, OUOE+Ol 

1,054E402 

l,243t.*02 

1,4076+02 

1,72iE+0? 

1.947E+02 

2,39U*o2 



2,076F*02 

3 1i;6E40? 

3,7<ibE*02 

a, lbSfe402 







U 

7,)i«Eq»i 

7^M?F4fll 

7.377F*01 

7,9U9E*01 

9,0136+01 

1,0U6E*02 

1,2306+02 

1, 4026+02 

1,6696+02 

1,9156+02 

2,1666+02 


2,OdljE402 

2^794F4ft2 

3,ln7E*0? 

3,42bE*02 

3.7S5E402 







5 

7,lt<IE*ai 

7 1 l2F40t 

7*3i4r4n| 

7,937E;01 

6.976E+01 

1,039^+02 

1,2126+0? 

1 ,4016+82 

1,6076+02 

1.8376+02 

2,0B6E*02 


2.346EtO? 

2 BOSF402 

2,8b2F40? 

3.126£*02 

3,3916+02 







6 

7,|lll£*(li 

7,1 ISF 4 OI 

7,3i2f*ol 

7*922fc*0J 

8,930t4ol 

1,0286*02 

1,1866*02 

1.357E+02 

1,5516*02 

1.7656+02 

l,99|t*02 


2,?l5ft0j 

a.OSflFtOa 

2 6B7r402 

2,698E*02 

3,1276402 







7 

7,|l‘lE*0l 

7,11SF*0J 

7,3ft9t*oi 

7,90SE*.01 

6,8786401 

l,0l4E+02 

1,1566*02 

1.3186*02 

1,5026+02 

1.699E+02 

1,9006+02 


2,07EE*07 

?,29BF*02 

?,5(>3E402 

2.707E+02 

2,9076+02 







6 

7,1I4E*0| 

7,il4F+0l 

7,5n7E+01 

7,887E*01 

6,6|9E40I 

9,96b£+ni 

1,1306*02 

1.2856+02 

1,4576+02 

1,6366+02 

1,6116+02 


1 ,‘ia'lE + Oj 

2 170F*r>2 

2,5slE402 

2.526F+02 

2,691£+o2 







9 

7,|laE+0t 

7*11«F+0! 

7 Jl3£+0t 

7*865£401 

6,7356401 

9,797E+0l 

l,t08E+02 

1,2556+02 

1.4136+0? 

1.5736+02 

1,7296+02 



2,04BF402 

2 2H3E402 

2,347EIo2 

2.479E402 







10 

7,1li|E*0} 

7, n4F*ftJ 

T*?49E401 

7,6386*01 

B,6«764oI 

9,65aE+0i 

1,0896*02 

1,2276+02 

t,3T|6*02 

1.512E+02 

1,6536+02 


l,79OE*0? 

l^RSlFtOP 

2,0s9£*02 

2, 174jeIo2 

2,2796402 







11 

7,)UE*0l 

7 1 1SE401 

7,2e3£*01 

7,607E;oi 

6,563E4oI 

9,5«1E401 

1,071E*02 

1,200E+02 

1,3296+02 

1.455E+02 

1,5816+01 


l,70ilEtCl7 

l^Rl9F*ft2 

1 ,flj3€*03 

2,0146402 

2,0976403 







u 

7.)lilE*0j 

7*ll2r401 

7 2A8E401 

7 774E;o1 

6,0956401 

9,4436+01 

1|056E*02 

|,174E+02 

l,290E+02 

1,4036+02 

I, 5156+02 



1 717F*0? 

1 BolCfO? 

|,875t*02 

1,9396402 







iS 

7.1UE*»t 

7.n0F401 

7,2|8E401 

7,74iE*01 

6,4416401 

9,3S9E+0l 

1,0416*02 

1.150E+02 

1,2546*02 

1,3576+02 

t«456i+02 

to 

l,54HE402 

1^B29F*02 

l,6«7EfU2 

l.TSBExOa 

1,8096+02 







7,|laE*CI1 

T,tftOF* 0 | 

7.270E401 

7,713E;o| 

a,S99E+ot 

9,268E+01 

1,0296*02 

l,129E+02 

1,2246+02 

1,3186+02 

1,4066+02 

15 

l,<j9fcE*0f 

J^B'jSFtOa 

1 6i3E*02 

1,664£*02 

1,7096+02 



7.llllE«0t 

).a>SE«07 

7, 1 07F*fll 
1 ,aOBF*02 

7,2b?E4.01 
1 ,Sm9E*02 

7.692F40I 

l.593£;02 

6,3696+01 

1,6306+0? 

9,2306+01 

1*0186*02 

1*1116+02 

1 20OE+02 

l,26bE+02 

1.36T£+02 

U 

7. 1 laEiOi 
IjUOJEtOj 

7^ToSF*(M 

l^fl«;7F»fl2 

7 ?«|bE401 
1 Srt3E*02 

7,o76t*0| 
1 SbSE^Oa 

6,3436+01 
1, 5776+02 

9,1676+01 

1 , 010 E *02 

1, 0986+02 

1,1826+02 

1,2636+02 

1,3386+02 

17 

7.|tHE«0l 

7^if|SF*fli 

7 Pii?F^01 

7,669E*0l 

6,3306+qI 

9,16lE+fil 

1,0056*0? 

1 ,0906+02 

1.1716+0? 

1,2496+02 

1,3206+02 

18 

l.38;E*0? 

KUS3F402 

1 496 E 4 O 2 

1,514E4.02 

!,5u6i; + o2 





7,tl0£t0j 

1,37EE»0? 

7,ififlF*01 
1 a?5F*02 

7,?qiF4ei 
1 4a?E*02 

r,6fc7E*.0l 

I«50SE402 

8,3256+01 

1.5356+02 

9,1S2E+0| 

1,0036+02 

1,0876+02 

1 1676+02 

1,2446+02 

1,3156+02 


Figure X2 - Continued 



THREE'OlHENSlONtL RllINT-HnDV fLCu (CFD.CnDEri) 

initial CONPITION^ HACH no. to' 00000 P.1«|F. 7.9T7N1E*02 RHO.INF* I , OJOOOf .06 CAHHA.INFA I.OOOOO 

ALPHA. n’ooOOOE.OlthFGI AINF* J.JPTPSEtOU OiNF. J.29705E.05 
NITEN. 0 NENO. ?(l nGAS. .) IL. tP2 JL* I2 AL. 10 JSHK. 10 


spo. I'.ooplE.oi r?« N oOnOE-oi sj« o, zfocnl. o.ptWEaoe 


NITER* 1 EIGENVaIUF TNFoI (rEC EIG.HAXrli T>K) 

-< I.SSoBE.OS.I* IN, 2. IrK* t)( O.OPISE-O*, I. t,J. 9,K. 0>( «,S«SI)f-il(r, I« 1,J. 

the MAx> a,l966E*0N CFL» S.tplJF.Ol OtC» a,5aeOE"00 0R« E.OSlItEAOO 01. 5.6W5E-06 
X(I,IL,JSHK). 9,6O0NE*01 X(Kl ,IL,JSHK)..J,^S0EE*02 E ( | , I , J*HK)».1 ,B091E*01 
P(1|H). i.^^ANETOj P(KL-i.K). j' 511PE.03 p( 1 , 1 )> 1 .OlfllEAOS H»lT7 


iiK. tr> 


HAX ENTHALPV ESBOB . •A.55tE-02 IC. IS JC. ^ Rc« 1 

HP, OF £BROB»-< II>i-0 * . 22«S - — — 

NO, OF errors < 1,1! * ■ 1J2a 

-NO, OF error* . 0,l-« . 12 o ----- — -- - - 


-llIUR. — 2._ElCeN»*l yp TNFfll-<REO-61C-KAX,-l,j,4(7 

I l,R001E.05,t. 1S,J. S,K* ?)( R.BIORE.OO-I. 1,J« ’,KR IK R,57JRE-0B,I» l,J« 1,K» 
THE MAX. J,!b<.SE*OS- CFL« 5,Ami 7TE«<H 0TC» «,ST3RE-0B- -OR*- *,«S6*£*«0 0T» StBiRlE-Oi 

X(1,IL,J8HK). 0 BlOREtOl X(KL,U,JSHK)».S,^305E*02 Z( I , I , JSHX)«.1 ,BSo 9E*01 
-1>14,XJ. S^ilpSE.OX— e<XL»ifM).-l 6B22E.03 p( 1, | R,7AsSE*Pa H»14T - 


MAX. EHTMAUtt-IRROB -P— • 

HO, OF errors a 10,0 X 
NO, oF ERRORS <-1,5 *-* 
NO, OF errors < O.I * » 


,6SiE.«2 Kiu. 

22oS 

*2. - -- 
SO 


-M — JC. »._Ke.__|. 


TTtER. 3 eigenvalue INFoI tREC £IO»HAX,l,j,K) 

_( |,3S23r.05,E« LW. -*,i(Jt_^)4-«.rHM£.W'l4. — 1W» — S,K*- tW. I-rX. — M 

THE MAX. 2,I731E«0S CFL* S.BBORE.Ol BTC* g.60lTE*OB 0R» 2,0530E«00 OT. 5.B602E.06 

.Xt4,IL,JSHKJ..0,5*47F4.Ol J1<K4 ,U^JSHI0..3 t. 210E*02 -Z( l, l,JSHKJ..l .SRBIEAO t - - 

Ptl.H). T,6E7RE*03 P<KL*TiH). 1 tO|7E*03 p(l,l). R,335lE+0« M»177 


HAX enthalpy error > •e.RRlE.oa 
— ERROR* a-?^,0-*-.-22oS 


|R JC. 


5 


KO, or ERRORS 
^ -Na^DF>eRRoAS. 


< 1.0 X • Qli^ 


■■^JITERm g— — - 

t I |}19SE*0S«1* 15|J* 8 »K*i P( 9*5859E*08.T« iH a*625SE«09»U l#J» 1#K» 55 

THIg. HAvt->, 1 a.i9£aq^ cgL-»--5r^i^TK*^ 0U*-Sw*^iUC»04 — - 

X(l,IL,JBHIt). 9,S73aE*OI X(Rl,:L,JSHK).. 3,R113E402 E(l,l,J0HAl«.l,84t8E«01 

PtlaHJ. ^,928tE*03_jsIXL»I,H}|,-t-5J46EA(L3 r( 1,U. B^9S37E*0«. M»l7J- 


HJX.ENIH*Lpa_ERBOR-J~«i,04olL.W -IC» — 1R--JC.- . A—Kc* 2 

NO, OF ERRORS < Ip'o X . 2263 

MOr. O E -ERHOM .. 1,5 K 3R» 

NO, OF errors < 0,1 « ■ 2S 


MITE*. S EiGEMVaLUE INFpl (REC EIG.HAX, I, J,K) 

U J-2a«2E.0Sa*-3iW* -SH*.— +iL-9.4325E.04^— 2,J»-9.*. ITU RaARROl.O*, 3* -Atjt -1. 
the HAX» 2.1S06E.05 CFL« 5,B(i 14E»01 BTC a.OORRE.OB OR* 2,035IE.OO BT. 3.4I0SE.04 

JUL,IL,JSHKJ._R,S444E+01- JL<JtL-,H.,JS**RK>«a.iOA‘C»2- -ltl..iW*MMC.A.,*3A8E.04 

Ptl,H)« l,2473E*0a P(KL.i,Hl. 1.544TE*03 4(1,11. 8,*29SE»00 H.jTT 


|T> 


NITER.S96 eigenvalue INFOI (pEC EIG-HAX, 1 , J,K1 

I I,216BE>0S,]. In j. 9,K< IK 9 S079£-1)4,I. 1,J. l,xa 2K 3.7319E»04,I. l,J. 2,K. 1T1 

the HAX* 2.4742E+05 CFl” 9 OxlOf-Ol BTC 3,7339E-04 OR" I.OSIOEaOO BT* 2 B»70E«04 
X(1,1L,J8HK)» 9.afc32E+01 X(Hl , IL, JSHX)a.3.a7S2£*02 l(l , l,JSHS)..9,aS8BE*0l) 

P(1,H). 2,860TE+0u P(FL*1,H'« 3.4(i 09E*00 4(1, 1)« R.SlSBETOd H.177 

HAX enthalpy error . O.RU7E-(i2 IC 1 JC. 2 Kc« II 
NO, OF ERRORS t 1o' 0 X . 2?oS 
HO. OF ERRORS <103° 2Z3o 
NO. OF ERRORS < 0.1 X P I03fc 


Figure 32 - continued 



HITcao<i<>T tirpNi/iiiif iNii'i (Ur ftc-nx,],j,K) 

( 1 .J 109 E- 5 S, Is (c, Js (,tc ,11 II SuRSt-lls, Is liJs l,Ks g)( J 7J3ilE-06, Is t<Js ZiXs i7> 

tHF MAx= 2.67SSEK1S CFl s 1 1 ,H 0 F.ai I.ILS 3,7SJSI-III1 OB' I OSlEE+liO Uts 2 9S77E.HS 
Xdrtl .J7HK)s B.ihnijf *r>l < ( k 1 , ) I , J9HK ) =-3 , a 750t I(l,l.JSI'k)s,9,ao0be»00 

PtliHJs 0(KI.-1,I1; J nkijinoa pf),j)s 9 b2SBt*04 "S 177 

MAX fUjH*l PV IRpnii s a ilpfF-n? ITS > JCs ? Kca || 

SO, (IF FBkflHS < 1(1 n S ' JAoB 

SO, nF FOk(iPS < ) 0 » s 'P ?7 

SO, OF EBFOBS < (1 I X s lOSx 

I'ITERbsof Erocsviiiit iHFti (ptr eUi-»a*.IiJiK) 

( t.aiRBF-Ob.I' IR Js 6,SS |)( a 5tO?E-OB,Is I ,Je l,K' 2)1 1.7}j7£.06, Is l,Js 2,Ks )7) 

THE H**s 2,(,7B1E*0R CH' 9,0 1 »RF-fll OTf' 3,71371-00 0H= I,(l5l6e900 OTs a.SOBfiE-OO 
,JSHK)s 9 lIsoaFttll Tlkl ,n ,JSHMS-3,a757Et02 71 1 , 1 , JSHK )*««, ooaOFiOO 

P(I|H)S ;, 629317(14 P(KL'I,I)' 3 o6oaE*Ou p(l,l)' 9,51966469 S'l77 

MAX FNTHAl.PT EBROR s 0 3S3E-02 ICn I JCs 2 X£s t| 

SO, oF ERRORS < )0 6 * s 2295 

NO OF ERRORS < I 0 « s 223? 

SO, oF ERRORS < 0 I X s 1999 

NITFRs>; 99 EIGFNTiliic II Foi (sEC E lO-MAX, I , J,K) 

( I.2I97E-05.IS l 9 ,js b,K= DC 9,5tI6E-0(i,Is l,J' 1,S' 2)( 3,75396-64, Is l,Js 2,Ks I7> 

THE MAX' 2.678264(15 CFl' 9.6,366-01 OTCc 3.73396-64 DH» 1.05266466 OTs 2.90006-04 
X(I,IL.J3HK}i 9,9S?SE4D1 X<Ki ,u,JSHX)'-3, 97516402 Z ( | , 1 , JSHK)S.9,64S3E406 
P(l,M)s 2,81096*09 PCX6-l,li»E 5 44026*69 PC1,D> 9,50596*09 HSI77 

MAX fSTHALPT ERROR s 9.37hL.(l? ICs 1 JCs 2 KCa )1 

SO, OF FRRORS < lo'o X s 229R 
SO. OF ERRORS < 1,6 X ' 2255 

NO. OF ERRORS < 0 1 X s 1 on,, 

MITERsaOO EIGENVamIE ISFOc (pEr EIC-MAX,!,J,K) 

( 1 ,22a5E-6S. I' 15 . Js SiKP DC 9 5t39F-(16,Is l,Js |,Ks 2>C 3.7339E-04.1s l,Js 2.Ks |7) 

TME MAX' 2.6782E46S CFE' 9 OoSOE-Ol 0IC» 3.7339E-0P DR* I.0524E400 OTs 2.9619E-64 
X(l,IL.J 5 HS)s 9 . 9987 F 401 XCK| , IL.JSHKJs.l, (,7516462 2 ( 1 , 1 , JSMMs.9,973IE400 

PCl*H)s 2,79726460 PCKA-l.Hls 5.95996409 PCI,D* 9.5136E409 ,1*177 


Figure 12 - Continued 



0» 

0> 


Kt I fMThAlPY kkROeS JaCOlUMK) 

. I ? X M •» b ' B 

1 ?' 6*>7F*(J?-6,X67F-rt5«l «3‘iF-f'3 4,V0?E-0* l^iOOF-Oa-i 036t-05-i 

2 7^«6eE-ft3»P ?S6F-(vX-S -05 I,J7»t*03^2 588f-«5 « 

3 3,2??t-43-3^n39F-&3 1 , 27?F«03-3* 79Vfc.(>6-2,<?'iSf -04 Ub2PP-(*« t,ltlit-03 1,7T0E*0X„J 

4 3 69tE-03«t> 962F-04 I <i7iF«03 2 032^*03 5* ] «3P »04>i ^SbU*03 \ 499b-04 |,b0ie»03*9 

5 •9^607F*&4w2,7a?E*&X*>l a«xF -r 3-« *eoat-0i» U63l£"0i» J ,a5bE»Oi-4,iM4E*03»2 275£^03»9 

6 •7.3^7f-'i4^l fl?2F-C3-2 947F*0j«U3'?8P -fi3 H.StbE-'Oa t^^SbE-OS^S 767E-03-I 
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a 

01 


>. V*BUSlfS (P0l,iF-»«rLr, Tl.tT^ (]>) (SOLE, Pt‘t (K)» RAOUL OIStANCef tj)) 


W vAR Ro I ([oROkI ltCCLU**N) 

1?aOSo7S9 1»11 

1 -1 7MB«oa a kOOBaOS ? %<>7 F*ao l,/?3t*Cu ?.()20t*0a 3 3Aof*ni ij,OiaE*AA il hRortOA S.6'l7f*0U O.JOOEAOi 3,?OflFA0S 

2 .<7jaE*(i3 1 ?03r»00 c (lAkEtoa ? A'3t40u J.^OOEiOo J.B76EtOO a,7nlE*tio 5,02»E*0« 6,7(j‘ir40« 3,J0eF»05 ),?<)8E40S 

3 K1?7E*iia J ^OfiEkfO 2.0hAE*(io R.aRRE'Ou O.IIaFaOi o,73oF-KIO S,JBo£*00 o^CROEaOA 0 7?0E»0o 3,JR6E+06 J.iOSEHS 

R a^BBOEkOU 3 (i^oPaO*! 1) “ 775E*l'i] ii,A53f4<‘0 k,07(,E*04 6,7RUEA<'0 7 307 Ea 0« 6 ??8E»Ca 3,SRI)E4I15 3,2»eE*l)5 

5 R,iMtEtou 5 ?07E+C'J 5 B'loriOA J,o65t«Co 7 ,bo 7E40« e.SOREAOO fl,flflaE«0« R,7JSf»(10 3.2RBf-»0B i,ZO8t40B 

* o 6J?C*<tR 7 ??3F»nu 7 77jiE<0R " 320F*Cia H.BBiFtRO o.357£tOa l,Ol5E«oa 1 007E403 1,I7<>F>0S 3,20 BEa 05 3,2RliE*(i4 

7 8’7*i7etna 0 jkhFACa R oOrFaDO IjOlRF^Oa l,r,76t»0i |,i27E*05 l.lBfE*OS 1 ZElEtOi l,33RFtl>'i 3,2>ltlE*0S 3,29eE«06 

» l,12«E*n5 1 ikoFtOS 1 lB7F40b ),??7E*03 l,2tsF»0i> | 3SBE40S l,3a»E*05 I.USBEaOS l.aSREABS S.aRBEAOS 3,2oaE4D5 

9 1 372£*oB I 39tF*l)S 1 3°(.F ,<15 l,o?4Ei03 1,939 Fa 05 |,SI2E«0B l.BJOEfOS I.oOOEaOS l,b93F*05 3,29BE405 3,298C*0S 

10 !>7ifto5 1 A18F403 1 oSRFtra 1,61(>F«1>8 i, 631F+05 | o54E*0b 1,713E + 0S t.7a?E»0S l,8f,7F«0S 3,29aE405 3,?98 Ea05 

11 l,997Eto5 1,9<I7F»05 l,6b<)E*0S 1.6UE*riB lieSZE+Ob t,787E«05 l,648E+05 l.8Il£*04 2,0a3E4CS A.SRaEtOS S.ZOBEtoS 

11 ?'S06E*0b 2,il9F405 ?,07kE*oS 1,903F»05 2,003E+08 l.R3«E*0b a.OOlE-tOS |.925E»05 2.130E405 3,Z98F»0S 3,29aE408 

13 2'3(lcFto5 2 ,?ii 5E405 ? SOtEtoS 2,I?S>E.08 2,1?9E4«5 2.O68E405 2,10oE405 2,oaoE4(l5 2,S25F40S 3,29BE405 3,298E405 

:« 2,357E405 2;i19E405 F,29iiF 4I)5 2.227E408 2,227E40S. 2,|60E40i> 2,2llie405 2, 119E«05 2.291F405 3,?98E«05 3.298F405 

IS 2'^a|«E40S 2 3BBF405 2 SUoFfOS 2 J86E405 2,2B9E40S 2.226E40S ?,277E40S Z.lRaE+OS 2 330F40S 5,29BE40S 3,298E40S 

FNr ViRUBIES (POLAR. AkrU, THETA (Dl BERlOlOHAL ANBUE, PHI (KJl RACIAL DISTANCE. (31) 

FhO v*fl Ks ) (1=R0H1 JcCOLUHN) 

123RSS789 t0ll 

1 3'9l6f41(l 0 IOoEaIO 9 OSBFfto U,{i76t4l0 a,0RSE4l(l 9,032E4l(l 9,9lSEt!0 R.OOREaIO 3.970F410 1,9R2E409 1.992E409 

I 0,0238410 0 (iBIEaIO 0 OlllFflO 9 03(8410 5,9988410 3,99(E4lo J,9S1F4I0 S.oSOEaU 3.8SOE4I0 1,902E40» 1,9028409 

3 3,9fl5tt|0 3 9P1F410 3 937F4IO 3,901E41f, 3,9(18410 3,8898410 S,(>S7E4lO 3 8208410 3,805F4((I 1,992E409 1,9028409 

0 3,88«8410 S.BSOEaIO 3 8587418 3,810t4lo 3,78oE4lO 3.7518410 3,70OE4lO 3 8S9EflO 3,8138418 1,902E489 1,9028409 

5 3 78AE410 3.735F4IO 3.7018410 3,4.81(410 5,8loE4lO 3,5Bo£4tO 3,5388410 3.O03E41O 3,0?lF4lo 1.902E409 1.902E409 

6 3'5BSE4I0 3,557E4lO 3.558E410 3;o92E4lO J,(l38E4tO 3,3918410 3,320E4I0 3,2928410 3,178F410 1,9028409 1,9028409 

7 3,OUE410 3 37BF410 3 30lE4l0 3.300E410 3,20o84l0 3,1838410 3.100E410 3,080£4lO 2.92IF410 1.902E409 1.902E409 

8 3,1B|E410 3 159F4I0 3 l?3(4lO 3. 0781410 5,038E4l0 2 941E410 2,9o9E410 2 B11E410 2.700E4to l,902F409 1.902E+09 

9 2'905E410 2,9718410 2 95of4lO 2,8918410 2.B538410 2.765F410 2,7228410 2,5M£4tO 2.55BE4I0 1,9028409 1.902E409 

10 2,7i«E 410 2 705F410 2.7UpF4lO 2. 7088410 2,65o84lO ?,59lEtlO 2,Sl5E4lO 2.0528419 2.3118410 1,9028409 1,9028409 

M 2,0398410 2,B08F410 2 5228410 2.501E410 2,o7aE4lO 2,0088410 2,530F410 2,3O7E4l0 2 107F410 l,902C4O9 1.9028409 

12 2,252E410 2,2918410 2,3l9FtlO 2,379Etlo 2,3198410 2 3278410 2,1888410 2,211E4lO 1.908E410 1,9028409 1,9028409 

13 2'iioE4!0 2 toiEtU 2 17n84lO 2,236E4tO 2.l9lE4lO 2,226EfU 2,078E4lO E.OBSEtlO l.aot£4lo 1,9028409 1.902E409 

10 2 Toio 84|0 2,0iO0F4lO 2,0098410 2,n2E4to 2.080E410 2,110E4l0 1,9698410 t,973E4lO 1.693F4I0 1.902E409 1.902E409 

IS 1 931E4IO I QBoE.lO 1 970E4IO 2,038(410 1,9838410 2 032E410 1,89BE410 I,902E419 1.8208410 1.902E409 1,902E409 

PRE9 VARIABLES (POL>F-A'<rl t, THFTA (Uj HEHIDIONAL ANCLE, PHI (K)l RADIAL DISTANCE, (J)) 

PRES vAr Kc I (laRObI JcCdIUHN) 

, 1 23056789 10 11 

1 9,5028400 1 (1108405 I 03bE* 05 1.012E405 l,OO7E40S 9,9018400 9, 6138400 9.662E400 9.0(10E400 7,9788402 7,976E402 

2 9’T35E40U 1 007E40B 9 97(iE,oa 9.B77F40O 9,8018400 917338400 9.079E400 9 5908400 8.627F400 7.978E402 7.978E402 

3 9,329E40U 9 o» 7E40O 9,079E400 9.010E400 9,38SE4C0 9,2928400 9,2278400 9,02OE.0O 9,061F400 7,9788402 7,976E4C2 

0 8,60S€404 8,707E400 6,70kF400 E 7|98t00 8.709E400 8.660E400 8,6058400 8,46BE400 8,0178400 7.978E402 7,9788402 

5 7'7898400 7,890F400 7 OSof.OO 7.9O|E400 8,0)2E40o 8.027E40O 6,002E40O 8 OIIE4OO 7,9SBE40O 7,9788402 7,9788402 

6 6,7B6E40a 6>22Et00 7 03NF4OO 7 1208400 7,l7sE40O 7,3068400 7,1896400 7,086E400 7.0658400 7.976E402 7,9788402 

7 5,6638400 5,0368400 5 9858400 6.130E40O 6.2“9840o 6.0218400 6.003E400 6^858400 6,5078400 7,9(86402 7,9768402 

8 0,5838400 0 000F4OU 0 9978,00 5,I71E40o 5,0008400 S,o9o6400 5,7988400 5 7628400 6,1658*00 7,9788*02 7,9788*02 

9 5j6t3E*flO 5 S6Sf*00 0 10oF*«J 0 312(400 0,6568*00 0,6696*00 5,2078*00 O.SR3E*00 5.8178*00 7,9788*02 7,9788*02 

10 2 7196*00 3,0S3F*OO 3,3566*00 3.63i(*Oo l,95o6*0o 0,1606*00 0,0318*00 O,7B0(*0» 0,05Sf*00 7,9788*02 7,9788*02 

11 1 777E400 2.159F400 2 00364(10 2,8878*00 3,1338*00 3,5968*00 3,7168*00 0 0938*00 3,7638*00 7,978£*02 7,9786*02 

12 i;20(lE*00 1 32of400 1 6I(|F400 2,23|E*0o 2,051E*00 5,012E*00 3,1006*00 3 967E*00 3,5796400 7,9788402 7,9766+02 

13 9'879E*05 1,1668*00 1 3936400 1 7598400 1,9788400 2,5158400 2,6526400 3.0888*00 3,1266*00 7,9768*02 T,978E*02 

10 8 5718403 1 0018*00 1 1S6F*00 1,0578*00 1,6“9E*00 2,l?2£*0o 2,2066*00 3 0918*011 2,S76F*0o 7,9768*02 7,9788*02 

15 7 2708*03 8 B62E*03 1,0828*00 l,269E*0a i.o3lt*Oo j 8536*00 2,0268*00 2 7618*00 2,7908400 7,9788*02 /.978E402 
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M*Ch V4PT*Plp*r fPPM 4P-4»‘rLf , Uif T* (M» *>*ClE» f'Ml {k)j «AUIAL. DISTANCE, <J)1 

HACH vAH Ka 1 flsROw; JsCui"Mri) 

i?-^tiSb7e9jon 

1 ?'572E-01 7.7t0f-0S I 7kfF-01 1 ?61f-Ol 1 d<>3E-0l a.JU'iF-ll) S 2?S£-(U 1.I1Q7F-OI l,0Ci(ie*«l l,OOi)fc*OI 

? l,77<if.(H l.PtOF.m F.SSsl'M E.fcSfct-Ol 3,J(>3E.n| i,u5it-OI 3,9j,>t-0| a I33E.III 4,l6if»(ll l.oOOttOl l.OCOEtOl 

3 3^?l5F»i't 3 7?lf*01 «,l*i7F*01 9,?02E*^l a 69jf-01 u 926E-0) E,30<>F-01 5 6R6E-01 6,0EUF*'>| l^OOOF+Ol I,000E+C1 

9 5^079C-ftl 5 jj^oFaO) 9»69^r,.0] 3 P/iwCf-Pl e 207E-P1 6 996F-PI 7*J09t-01 7 ^10E«01 l*000E+01 l,OOOEfnl 

5 6l599F-»l E,fll)lF>0l 7 lipE.Pl ;,i7ir-El 7.9e3E-<M B 33OE-0I 9 12BE.01 9 9olf-pi l,0O0E«O| l,090t*El 

4 e'319f»0l B,f,F7E>l)l “.BUhE.'il 9,J39E-Ol 9 71lF*01 1,0 B\jF 409 1,059F«FI) l,|49E«0t 1,19EF*B9 t,000E40| 

7 9,6?St-Pl t ot5F»0(i 1 0'lf.F + iP 1 C87E*Pp I.I3BE + 00 j iVjttOO |,237t*00 l,?90Ftl)0 I.372FEBI1 l,OOOF*OI l.OBOEfDl 

B )'l7?En(E ),t93F»PE I 2FPF4-I0 l,e6itE*0li l,304C*FO 1 37bE*00 1,413F«00 l,500E*tl0 I 5?4E»0S 1, 0004*01 (.OOCEtOl 

9 I 3iut*00 1 3B9F«00 l,3FlE*00 1, 9151*00 l.PSaEtOO 1 S33E*00 1,595F*00 l,710F*00 l,7«2f*OI> t,000F*O| t.OOOE*OI 

10 1,6B6F*00 l,SblE*00 I 54FF*00 1,589E*0(. 1,427€*00 |,48 iE* 00 |,787E*00 ),833E*I)0 I.994E*0(I 1,OOOE*01 l,OOOE*Ol 

11 I 897E*00 l,7F9C*00 1,77;F«00 |,79Bf*0l> l,799E*00 1 808E*0l) l,93df*00 I,9O2E*0O i,20U*00 l.OOOEtOl 1(OOOE*01 

1? 2'ou8f*00 I 990F*00 1,96<iF,flO l.VOBEtOo l,9MEtOO 1.990E*00 F,045E*00 2 044E*00 2,356f*0» 1,OOOF*I11 l.OOOE+Ol 

13 2'tSBE*00 ?,13SF*00 ? 09(iF*00 2 OOSEtOO 2,0844*01) 2,U59E*00 2,181E*00 2 183t*00 2,509E*00 1,OOOE*01 1,OOOE*OI 

1“ 2,2594*00 2,?32F*no 2 2!iiE+00 2,165F*00 2,205f*00 2,l8lE*00 2,illF*00 2,3104*00 2,6304*00 l,000f*01 1,0004*01 

15 2 395r*00 2 ?Of,F*no 2.293F*«0 2,238F*C0 2,289E*00 J 257E*00 2.399E*00 2,397E*00 2.720F+O0 l,000£*01 l,OOOE*Ol 

PHO VOFIABIES rPOLAF-Af'cLE , 7H4*» (I)( HEBIUIOBAI iNSlF, PHI <K)l BAOIAL UlSfAHCt, (J)l 

BHO VAR Do 2 noRO)«l <:C0UIHM 

IJ39547B9 10 11 

1 6'197E-06 S.fl19f-0b b bOoE-Ob 6 150F-Ob 4,269E-06 b.lblE'Ob 6,122E»06 b,0O9F.06 5,925F>0b 1,027E-Ob l,027E-04 

2 4,097E-04 6,09UE,06 h.lOpE.fib b l?7E,0b 6.1124-06 6,09eE-Ob 5,9694-06 6,O63E-06 S,58IE*0b 1,027E-06 l,027E-0b 

3 s' 9936-06 5,963F-0b I 00?F-0o 6,0064-06 5,9924-06 5,9654-06 5,963E-06 5,881E-06 5,9534-06 l,027F-06 1,0274-06 

9 s'6I2E-06 5,6BoP- 06 5 72BF.06 5,7734-06 5 7924-06 5.826E-06 5,8l9E-06 5,899E«06 S,760F-06 l,027E-06 1.027F-0b 

5 S 2I5E-06 5 30bf-06 5 38|E-06 5,965E-06 S,530F-06 5,629E-06 5,629E-06 5,B32£-06 5,5624*06 1,0274-06 1,0274-06 

6 «'690E-06 a,fllo£,0b u.93iF-06 5,055F-06 5,l9lE.06 5,327F-06 5,«37E-0b 5,6094.06 5,5774-06 1,0274-06 1,0274*06 

7 «'^125£*06 9,2874-06 o 95)E-06 9.615E-06 »,B22E-06 9,9994-06 S,307E-06 5,2354.06 5,B52E-06 1,0274-06 l,027E-06 

8 3,593E-06 3 759F-06 3 96aE-06 9,198E-04 9,9164-06 9.6674-06 9,990E-06 5.230E-06 5.9224*06 I, 0274-06 1,0274*06 

9 2,827f-06 3 o93Er06 3,32iF.o6 3.6634-06 5,8784-06 9,3164-06 9.962E-06 5.12OF-06 5. 0354-06 1,0274-06 1,0274*06 

10 2 1624-06 2,9994-06 2 bBpF-Ob 3,0584-06 3.298E-06 3.780E-06 9,0204.06 9,7104.06 5,3264*06 l,027E-06 1,0274*06 

11 11706E-06 1,0994-06 2,2974-06 2,5234-06 2,BS5E-06 3,1824-06 3.7914-06 9,1334-06 5.6394-06 1,0274-06 1,0274*06 

12 1,596E*o6 1,7554.06 2.0194.06 2,2064.06 2,6054.06 2,8504*06 3,551E-06 3,8984-06 5,3954-06 1,0274-06 1,0274-06 

13 I’95|E*06 1,5954.06 1,8924-06 1,9894.06 2,3»84-0b 2 6994-06 3,3894*06 3 6824.06 9,6964-06 1,0274-06 1,0274-06 

19 1,3224-06 1,9984-06 1,7354.06 1,8634-06 2,2804-06 2,5534-06 3,2534-06 3.7504.06 9,336E-06 1,0274-06 1,0274-06 

15 |'ii2E.06 1,3914.06 1 5674-06 1,7914.06 2,1014-06 2,3894-06 3,0164-06 3 6134.06 9,9954-06 1,0274-06 1,0274-06 

U variables (PDLAP-ANrLB, THETA (U) HERIDI06AL ANCLE, PKI (Kl| RADIAL DI8TANCF, (J>) 


8 VAR K« 2 (IRRORI JoroLUHNl 

123956789 10 11 

1 l' 9234*03 8 9014*03 7.87(i4*03 7,9124*03 t,093E*09 1,0524*09 l,2l5E*09 I 2164*09 1 , 9724*09 0, 0, 

2 2,2924*09 2 955E*09 3 50(l4«09 3,2354*09 3,9764*09 3,3724*09 3,362E*09 3 9524*09 3 3324*09 0, 0, 

3 0,7614*09 9 665E*09 5,1034*99 5,0936*09 5,2974*09 S,3t0-*09 S,928E*0« 5 9874*09 5,o9Be*C9 0, 0, 

0 6,6?2E*09 6,d 9|F*09 6.8Bb4,09 6.9064*09 7,0964*09 7,192E*09 7,3794*09 7 5364*09 7 6324*09 0. 0, 

5 6, 13BF*09 fl 1684*09 P,279F+09 8,3554*09 8,9894*09 6,6594*09 ft,779E*09 9,190E*09 9,0024*09 0, 0, 

6 9,177E*09 9 |86F*09 9,26bE*o« 9,3664*09 9,5754*09 9,766F*09 9,8864*00 1,0|9E*05 1.0174*05 0, 0, 

7 9,B05E*09 9,7866*09 9.8674*09 1,0004*05 1,029E*05 l,053E*0b 1,0724*05 1,066E*05 1,1394*05 0, 0. 

8 9,9394*09 9 8954*09 1 0ObE*05 1 0464*05 1,0“74«05 1 065E*05 1,1006*05 I 1174*05 1.1594*05 0, 0, 

9 9,1834*09 9.9?tF*09 9 6124*09 9 950£*09 1,0184*05 1.066F+05 1,0764*05 1,1986*05 l,l|9f*0S «, 0, 

10 7,7826*09 e,962f*09 8.597E*09 9,1724*09 9,3314*09 | 0096*05 1,0074*05 1.1106*05 1. 1154*05 0, 0, 

11 6C7674+09 7,9894*09 7 57qF*09 6,1074*09 8,9226*09 9,1036*09 4,9994*09 9,9706+09 1,1164*05 0, 0, 

12 6'330E*09 b 7974*09 6 VS(l4*0» 7,3264 *09 7,7696*09 6,1236*09 9,V|84*CI9 9 9354*09 1,076F*95 0, 0, 

13 6,0794*09 6 1850*09 6 99|F*09 6,6584*09 7,2296*09 7,3296*09 8,9824*09 6,3954*09 9 8564*09 0, 0, 

19 5 7686*09 5 0894*89 6 ltoF,09 6 1266*09 6,7596*09 o, 7696*09 7,8936*09 7 8914*09 8,8664*09 0, 0, 
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VAfltAftlf P vrfTrP C( f fUt i/Jl <K>J W/PJAI MSrA»jft, #J>; 

¥9 I A*tnO»(jJMTF 1a^^^ll^•^) 

I ? X C *1 b 7 e q 10 J 1 J 2 

1 <• ■» p p ft, 5 0 y# o» 0* 

2 7 U3?e<oo 7 vi?r*0v 7 7 7hHt»n^ 7,(»7ut*00 7 '>85£»0<* ^.O^oG + flC 6»?0ot+00 ^*3l7r»00 » a2/f*00 d,!»56E + 0U 0,OH^E*O8 

3 J^^*7ftF40l \ S0PF401 \ S?sr»ll l,*^7?|+0l i.SyBF^OJ |,M5f40i 1 6«2t40j 1 bis5f*n l,7Wt40l l,735F*0l 

j ifl<*nun ? ?TnP401 ■» f'7,Uri a ,»,JiJ*»£4ftl 2 3ft2MfJ 2,ui9f-»01 ?,‘*‘it>E40l <! 4J92 F-»o| 2,5?9M0I 2,50bb4M a.bOJE^OJ 

5 2^?‘>2e401 ^ 2 ^ 3 ? iiijEfftJ 3,19J£4ej 3 3,2qjf*nj i,3yif+ft| S,39|E401 3,dftlE401 

A 3^5S5E + 0 J 3 i,7iH‘e»0l S.^CIEtOJ S,yt»Jf»Ol 3,^?&E*0l 3 98bf40| 9 Ou5F*fU «,j09F»ni u.UU^Ol «,23Jt4f}l 

7 «^179f*ni <J ?q?F*01 A %ast*0| 4 ^97fc*ii <j,ii7Tit*(‘i (J ?itJE*Ol a.OjOE+OI u tt 7MF4pj 4,p3ur4<ll u,Q07E401 M,q79E*0l 
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Figure 13.- Shock stand-off distance 
= 10, a'= 30°, ZFOCNL = 7, 
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Figure 19.- Shock stand-off distance, = 10, a = 30°, 

ZFOCNL = 5. 



Figure 20.- Magnitude of maximum energy error. 
M = 10, a = 30°, ZFOCNL = 5. 
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-o.nSEtSl -3.6076*01 -U156t*01 - .BOOF.«01 -■) .334E.t0 1--9..flfl5e*Ol -1.0Z7.£tflZ- . 

-D 295£*0i -3 64St+01 -I 166E+01 

•1 ?Mt±.OJ 7.3i'lk£JiLL-=Sj21U£i^U_i9^9a0t*OX.=a..il6H±QJ_ri.,9Mttill.-liO35l:»l>2. 

-6.374t*Oi -3 683b*01 -1.176E+01 
-6 IBiiEtal -3.721c»01 -1 ia7E*01 

-; .;36E*oi_- 6. >7Pbtai_=i* .^Ai;*aJLjSU.3Jat±ai -9.im97*oi - i.fla^^a^ -r i. oa zEfjiz — 
-O-SaAE+Ol -3.769E*U1 fl.a4f,F.Qi 9 36ap.>ni 1 .nn7c«n? i.fi94 c *nj -i.TAnF^n? 

-6.61.<iEt3L -3~197fc*of -l.7J6t-.Jl 

■5.5 2 3t*jl - 6.9U F-E01 -3 1036*01 - 9 0396*01 -9.730E101 .-JliJ?6l-*0S-jt.Jt71E*07 -L»U9E*0J - 

-6 691k-*0l"-3 6366*01 -1 2i9t*01 

t.fti.7P .nl -T.naiFtai -9. 3966»01 - 9 . 9106*0 1 r l .011 t*J)2 - 1 . fle . 9 i. t07 . j l.ia 7 £ t a2 — 
-6.779E*0l -3 6716*01 -1.23CE*03 

■6.B10F*m -7.?ljF«01 -6.192 6*01 -9..A1AiA0.1_^l-t.fllaL-tflZ -1 .0676 -707 -K l OTfinp - 1.1 5S6t0 6 — 
-6.6S3E+01 -3.9126*01 -1 2106*01 , 

5.913 -*01 -7.1726* 07 -9.1176401 -7.0376*07 -1 . nB6E±DZ^l,.-125£taZ--JL. llltt B2. 

-6.9"336+0l -a 9506*01 -1 2S1E*01 


2 


Flgute 25 - Continued 



STEP 

3C CTheU 

2 “C fjL+E> 

3(.0 IHrTA 

. 1*7 ^^7 

DEO 

SHOCK 







«.o02)t+C2 

i ( J ♦{ 2 

t 5<JAc02 

■» 

£ 3I9E402 


3«7«3&4>C2 

A 

4*4v4t+J2 

4*3i-7- + C2 


3C0Y 

»^ADIUS 







1 A7Sfc+02 

1 .-27E*CI2 

i«&lBE4U2 

1 7e«0EfO2 


2,406E*02 

2«S15E402 

2.577fc+02 

2«o06Et()2 

2»bl7L*02 

STEP 

31 DTHETA 

2"c&2t+C». 

DcVr TriETA 

dV*26a 

DbO 

SHOCK 

RADIUS 






2*024 EfC2 

2«CA3Ef02 

2 D99E4-02 

2dl94t^02 

2.328c »02 


3*8<.3E+C2 

4* 194E4C2 

4«59LE*02 

A«990E<»-02 

t> 334E402 

BODY 

RADIUS 






* 45VE+-02 

i A77k+C2 

lobilOl+Oa 

1 bZlE^OS 

1.754E402 


2.43bt^C? 

2 556t*02 

2o622ct02 

2.654C402 

2.66<tE40~2 







FLtM VARIABLES# NADV 

EOUALS 

31 










^.359e402 


Z»b87S«02 2,6a7E«02 3»069£+02 3.296E402 3«9^4£402 


2«6b6Et02 Z*bbbl^^Z 


OEHSITY 


12 

3 ><.9ie«»06 


13 14 

3> A66_E «_afe__ 3 >_^ ?06-D 
4«578E''07 3*113E-ar 


15 


16 


17 


4«5?8E''07 3*113£-()7 2*3SOE-*07 2»068E-07 2«33BE-07 

8.EE9E-07 6>120E-07 ‘■•(!>et-07 3.«72t-(;7 2.92<ie-37 2 917E-07 

_3 i«26VXTQL<'_ 

9t745E'^7 7.323E-07 5.a98E-07 «.2a2£-07 3,578E-07 3.279E-07 

_i 4.592 E-06 *.573E-06 9.520E-Q8 9.4‘i5E-06 9 377E-0 6 t.l90E-08 

1.2C0E-06 9 392E-C7 7.23OE-07' ? 67iE-07 9.677E-07 9,379E-a7 

A.7.62Er0O 9.8 81E-Llt. 9.53 2E-D7, .Aj.i2Jtl.-M.. 

1 95tE-06 1.J76E-C6 9,<,93E-07 7.7Q9E-87 6.317E-07 5.669E-97 

_ft_ iiSLij=j:£>_5:.fi.«.J:-r06 ±,325^3:6* -9.9092-86 9.722E-06 9.9a9E-o& 

1.789E706 l,999E-«6 1.298E-06 1.060E-06 9 195E-07 8.009E-07 

i««^t7E-n^_A.UaE=i>(»_.i .2t9E=n6_3,.Jll£=Ml 3.0186-llt- 2.721E-0 6 2.521E-0 6 

2.278E-06 1.980E-C6 1 7UE-J6 1.500E-06 1.323E-06 1.168e-06 


U25a6.-Jl6__Ii434£.-!!6_UIM£=aii._1.0MErOL4 

U3Z9E-06,_UZ23E=nt 

L:a&_.LtaJS£.HHL_JLiASl£z04_a.£45«E--96 

l- 04 .■uao.9Er.afi...i.6a4£=atL— 


9 .596C-06 

2 893E-06 2.t07E-06 2.316E-06 2.087E-06 3.889E-06 3.733E-06 

-3. LtXlSIreS LXgMrfit— S ■3-96.E-Q6 ■ 5 373E-06 9.3276-06 5 lBOE -06 5. Oajj-j; 

10 1^0276^0** 1^077?^***^ 3.J98E-S6 2.913E-U6 2.709E-08 2.590E-06 

"”l^027E-06 1.027E-C6 1. 0276-06 1.027E-06 i?ofn^6“^?027E-l)r 

'iJ—LJimzfiB • - ^ " 


-9.019E-fl6 UUlSEsOE 3.319F-06 3.ia?f-06 

■1,..027EtI 16 . l.ll27E-:l)6__1..02.7£r.l>A— I.ilZIE-at— 


l..£.22£:^86- -L.027Erl) 6 1 .0276-0^ 1.0 2 

1.627E-06 1.027E-06 1.027E-06 1.027L-06 l.O27E-0r'l^027E-»6 


Figure 25 - Concluded 



results DF IHTtkFKt FfTLEES BU'<r BDOr COOES 
AND STRFAMLIHE FETRIC FRGSRAN 

UNSTEADY bLURT-EOLY TAPE READ. 16 STATIONS ON IT. 

STEADY 6LUNT-B0CY TAPE READ. ‘t STATIONS ON IT. TOTAL MO. OF STATIONS ■ 19 


ALPHA * 30.0100 DEO 

RINF » .1993E-05 SLUGS/FT**3 

PINF ■= .lOfitFAOl L8/FT*»a 

KINF = 10.0000 

GAMMA » l.ACCOO 

Z/ZHAX-LOCATUINS REQUESTED 


1 

.0250 

2 

.ObOO 

3 

.075 0 

4 

.1000 

5 

.1250 

6 

.1500 

7 

.1750 

8 

.2000 

9 

.2250 

10 

.2500 

23 

.5750 

24 

.6000 

25 

.5250 

26 

.6500 

27 

.6750 

28 

.7000 

29 

.7250 

30 

.7500 

31 

.7750 

32 

.8000 

33 

.8250 

34 

.8500 

35 

.6750 

36 

.9000 

37 

.9250 

38 

.9600 

39 

.9750 

40 

1.0000 


ERROR IN BODY SURFACE DATA. STATION _.38<. DISREGARDED 


Figure 26 - Selected sample output from interface program 
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Z-STiTIGN *<0. 
4« 1 

2 Z= 





K 

PH 

Pi»£Sil«r 

UENilTY 

U-V£L 

V-VEl 

k-VEL 

1 

0. 

* S632E<C2 

X.X370E-05 

l,0b2«t+03 

9.4574E+02 

0. 

2 

X 1250t+01 

X.^'.’TE + OZ 

X.X329E-0J 

l*259U+03 

1.1037E+03 

2.93866+02 

e 

2.25CCE+il 

I E93XE+C2 

X,CO7jt-07 

l*3bl2E+03 

1.22JIS+03 

X.90726+C3 


3.3 7bOE+al 

X.7X33E+C2 

X.L897t-Ci 

1.O0356+03 

l,4382E+03 

X 7790E+03 

b 

-l.SOOUE+OX 

X.-596E+02 

«. 7329E-0O 

l967‘76+03 

l,52e5e+03 

2.'.938E+03 

6 

5.6250C+JX 

1.3691P+C2 

c*bZX5P-C6 

1.637bE403 

l*7077f+03 

2.89o3E+03 

7 

O.75CC6+01 

X.X687E+02 

7 V396E-06 

2.l802t*03 

2.06146+03 

3.299XE+03 

6 

Z.^ZJOE+OX 

9.733CE+CX 

6.9709E-06 

2*5218E+03 

2 4134fc+03 

3.o0o6E+03 

9 

O.OCCOE+31 

7.<-067E+CX 

, 9880 £-06 

2.9230.+03 

2.bl316+03 

3.7970E+C3 

10 

1.012iE+02 

t .2V87E+CX 

6.0727E-06 

3.3792=+03 

3.264LE+03 

3.7339E+03 

11 

1.125CL+3Z 

“i.9992F-*CX 

9.2666E-06 

3.854Cfc+03 

3 72686+03 

3.60X6E+C3 

12 

1.237>E-*02 

3 87l3k+OX 

3.EX76E-06 


4.1467E+03 

3.290X6+C3 

13 

1 350CE+02 

3.Xt29E+CX 

3.O799E-06 

4,713A£+03 

4,54036+03 

2.E659k+03 

14 

l.<i62iE*02 

2.tcX5EH0X 

2.6950E-06 

5.029^£+03 

4,8247=+03 

2.2857E+C3 

15 

X .5 7S0E+02 

2. X600P+CI 

2.3X0XE-0O 

5*3406E+03 

5.10316+03 

X.70O0E+03 

16 

X,6e75£+02 

2.Cb29EtOX 

2.2355E-U6 

5*4930t+03 

5.2328E+C3 

3.82356+02 

17 

l.tfOOOE+02 

2 Cfa29E+0X 

2.2355E-06 

5•53^«E■^03 

5.2650E+03 

X.7975E-05 

K 

Phi 

Cu2 

C62 

CNO 

CC 

CN 

1 

0. 

7.28I3E-03 

2.7377E-02 

0. 

0. 

0. 

2 

X.X25OH+0X 

7 2aX3F-03 

2.7J77E-02 

0. 

0. 

0. 

3 

2.25CCE+OX 

7.28X3E-C3 

2.7377E-02 

0« 

0* 

0. 

4 

3.37SCE+OX 

7.2OX3E-03 

2.7377=-02 

0. 

0. 

0. 

5 

^.5cooe+oi 

7.tbX3£-03 

2.7377E-02 

0, 

0. 

0. 

6 

5.62i0e*0X 

7 26I3E-G3 

Z.7377E-C2 

0. 

0. 

0. 

7 

6.7500E+0X 

7.26X3E-03 

2.7377E-02 

0. 

0. 

0. 

8 

7.8 7iO£+OX 

7.28X3E-03 

2.7377E-C2 

0. 

0. 

0. 

9 

9.0O00ct0X 

7.28I3E-03 

2.7377E-02 

0« 

0* 

0. 

10 

X.OX26E+02 

7.2BX3E-03 

2.7377E-02 

c. 

0. 

0. 

11 

1.X25CE+02 

7.26I3E-03 

2.7377E-02 

G. 

0. 

0. 

12 

X.237te*02 

7.28X3E-03 

2.7377E-02 

0. 

c. 

0. 

13 

X.350CE+02 

7.26X3E-C3 

2.7377E-02 

0« 

0. 

0. 

14 

X.962i.E-*02 

7.26X3E-03 

2.7377E-02 

0« 

0. 

0. 

15 

X.5750E+C2 

7.28X3E-83 

Z.7377E-02 

0. 

0. 

0. 

16 

X 687SiE+02 

7.Z8I3E-03 

2.7577E-02 

0« 

0. 

0. 

17 

X.8CC0E-tC2 

7.28X3E-C3 

2.7377E-02 

0« 

0. 

0. 


J» 2 


K 

PHI 

PRESStiPE 

DENSITY 

U-VEL 

V-VEL 

V.-V6L 

X 

0. 

1.9562E+02 

1.1985E-05 

1.3969E+03 

9.781BE+02 

0. 

2 

X.1250E+01 

1.9339E+02 

1.1399E-05 

1.5030E+03 

1.0739E+03 

6.7681E+02 

3 

2.25CCE+01 

1.B97CE+02 

1.1091E-05 

1.5600E+03 

1.1931E + 03 

1.3293E+03 

9 

3.3750E+O1 

1 72126+02 

I.0526C-05 

1.6670E+03 

1. 29036+03 

1.9366E+03 

5 

9.500CE+01 

1.5612E+02 

9.6555E-06 

1.8999E+03 

1.98596+03 

2.99996+03 

6 

5.62506+01 

1.3799E+C2 

9.0C62E-06 

2 0536E+03 

1.7609E+03 

2.9827E+03 

7 

6.7500E+01 

1.18796+02 

&.1529E-06 

2.39256+03 

2. 08296+03 

3.3962E+C3 

8 

7.8750E+0X 

9. 9379E+C1 

7.2159E-06 

2.6999E+03 

2.9833E+03 

3 6320E+03 

9 

9.0000E+01 

B.1130E+01 

6.28156-06 

3.10626+03 

2.9958E+03 

3. 77626+03 

1C 

l.O125t+0Z 

6.5159E+01 

5.9C38E-06 

3.5978E+03 

3.9Z79E+03 

3.7690E+03 

11 

1. 12506+02 

5.16796+01 

9.6ClB6-05 

3.9753E+03 

3.8860E+03 

3.6328E+03 

12 

1.2375E+02 

9.30916+01 

9.06796-06 

9 38556+03 

9.33696+03 

3.2992E+03 

13 

1.3500E+02 

3.67596+01 

3.5713E-06 

9.7593E+03 

9,69136+03 

2.8356E+03 

19 

1.9625E+02 

3.0722E+01 

3.2191E-06 

5.0959E+03 

9.98986+03 

2.29136+03 

15 

1.5750E+02 

2.9819E+01 

2 73596-06 

5.9319E+03 

5.2C606+03 

1.6539E+C3 


Pxgiire 26 — Continued 
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J» 11 


K 

PHJ 

0'<£SSLR^ 

CkNSltY 

O-V- L 

i-VFL 

W-VEL 

1 


i.duV3E+C2 

1 .l3^>6c-o5 

J.3i73r+03 

8* 22966+02 

0. 

2 

l*125ot+oi 

i.79iit4:2 

1*1 t-u 5 

3.32‘-oc + 03 

8*7A07t+C2 

3.b9v7=+02 

3 

2 250uc'4'0! 

1,73c1F4C2 

14x^37,-05 

3. 3711^+03 

1 .0AP5E+03 

1 • 72db6+Cj 


3.i7l,0\ *0i 

I*c525£+C<: 

1,L257£-^)5 

d *^o62E + 03 

1. 3133E+C3 

2.53AIE+C3 

5 

3 LCC^ 401 

1 ^^02E*C2 

I.IJ 30L-75 

3.baa<iE+03 

l*7049E+03 

3*222dc+C3 

6 

5,6cS0E+3i 

l,<.07br+02 

I.CvZZ'^-Od 

3.7o0ze+03 

2. 1720E+C3 

347713E+C3 

7 

6*7 j tC.- + 0 1 

l*2713Mo2 

i GtObf-35 

3 9500‘=*03 

2.7U66E+03 

A,18b3E+03 

8 


K126AE+C2 

l.C5A0£-wS 

9.17C6t+03 

3*31AAE+03 

9.9235E+C3 

9 

9*J(CCc-KTl 

9.t77?E*Cl 

1.01^Al-L5 

9.9<J76=+03 

3*9686E+03 

A4&013E+03 

10 

l.J12‘-r£4D2 

a 2o«3E40i 

<3. 7lA^g-06 

«.,797‘'E+03 

<V.6337E+03 

A. 3725E+G3 

11 

1.12*‘Oc + C2 

7.CA8A6+01 

9*2333C-06 

5.0705E +03 

5*2t09E+D3 

A*06b05+C3 

12 

I.>:7tt+J2 

b.e^ACE^Ci 

fi A»06E*Oo 

j.353fcf+03 

5*b056E+03 

3.589aE+Ci 

13 

1.36CCt+C2 

3.ir691f + Ci 

fc42<iA7t-C6 

S.391S++03 

6. 1312F+D3 

2, 9506!; +03 

1<> 

U^t25£ + 02 

J.*A653£ + Ci 

^.•233AE-ce 

S,2926++03 

o.9bi5t+03 

2.3033E+C3 

15 

14O75LS4-02 

i.392CE+ri 

** 4 i63Ac.-06 

b.2358=+03 

6*7eA3E+03 

L*6A33E+C3 

16 

l.fc57tE402 

i 43^716401 

F.^t37r-0t 

E.219',E+03 

O49A35E+03 

8.7752E+02 

17 


5,3S.3*E4C1 

b.3C33t-Cfc 

j.2208t+33 

6.9977E+D3 

1,7«77S-C5 


K 

PHI 

C02 

CN2 

CNQ 

cu 

CN 

1 

D. 

7.2813E-C3 

2,7i77E-02 

0* 

0* 

0» 

2 

i.i2iuc+01 

7.2«13E-C3 

2.73771.-02 

0« 

0. 

04 

3 

2.25UE+0I 

7.2613E-03 

2.7377E-02 

0* 

04 

0. 

A 

3.3750E+01 

7.2813E-C3 

2.7377E-02 

0* 

0« 

0« 

5 

9.SCCCE + 0J 

7.E813E-C’ 

2.7377E-C2 

c. 

c. 

3. 

6 

E.t)2t08 + 01 

7.2bl3‘^-CJ 

2.7377L-02 

0* 

04 

0. 

7 

64 ? 5CuC+j1 

7.2«13E-03 

2.73777-02 

0. 

0* 

0. 

6 

74«750E+J1 

7.2tl3t-03 

2.7377C-02 

C 4 

c. 

0* 

9 

94 OOOCF+OI 

7.2bl3E-03 

2 7377E-02 

c« 

0 

0. 

1C 

1 31256432 

7.2el3E-03 

2.7377E-0Z 

c. 

0* 

0. 

11 

14125CF4C2 

7.2bl3E-03 

2.7377E-02 

0« 

0* 

04 

12 

l*2375t+Qi 

7.2bl3E-03 

2.7377E-02 

0. 

0* 

0. 

13 

1*35006432 

7.2S13E-05 

2.7377E-02 

0. 

04 

04 

lA 

l*A62^t402 

7.2b13E-C3 

2.7377t-02 

04 

0* 

0. 

15 

1*57 iLc 4 j2 

7.2813E-C3 

2.7377E-02 

0* 

0* 

0* 

16 

l*687tt402 

7.2813E-03 

2.7377b-b2 

04 

04 

0. 

17 

l*dCaOE+32 

7.2U13C-C3 

2.7377E-02 

0. 

0. 

0« 


PB FROf' 6FC^3,Ke FROM INT FRPOL At IGN, THE RELATIVE ERROR* OEKIVAT IVFS H.R*T. Z ANO ?hi + RS 


K 

Eb 

PX 

RBO 

R8Z 

RBPHI 

RS 

1 

l.o913E+aO 

l.to27E+CC 

S4i093£-03 

d.73A5E-01 

0. 

2.1862E+00 

2 

1.691CE+bC 

1.6e38E+C0 

A.2370(.-03 

P.7705E-01 

-3.AA62E-03 

2.1919E+C0 

► 3 

l.'i9COE + OC 

1.6798E+CC 

64CA39E-O3 

E.8725E-01 

-5.3721E-03 

2.2065E+00 

A 

1.6HPbE + 00 

l.<.76r£ + 0Q 

7.3A65fc-03 

9.023SE-01 

-B.3A23E-03 

2.2268E+00 

5 

1.6bttL+00 

1.6721E+CC 

a.7172fc-0J 

9.2001E-01 

-9.07A6E-03 

2.265SE + 00 

6 

1.6Bt.C^+00 

1.6713E+00 

O.1712C-03 

9.3730E-01 

-b,A797E-03 

2. 30 706+00 

7 

1.6e3t£+d0 

l.cbAAE+00 

l*l3A9fc-02 

9.5150E-01 

-5.6681E-03 

2.3672E+C0 

a 

1.6S2SE+0C 

1.6613E+0C 

l4256UfC2 

9.6035E-01 

-3.927AE-03 

2.A335E+00 

9 

1.6b2CE+0C 

1.6S30E+00 

l472&3b-02 

9.6238E-01 

2.2393E-08 

2,51416+00 

10 

l.ol>27E+00 

l.tA31E+00 

2.2323E-02 

9.6198E-01 

6.A792E-03 

2.61206+00 

11 

l.boAbc +G0 

l.£A5AE+Cr 

2,32226-02 

9.fc0e5E-01 

1.2011E-02 

2.7339E+00 

12 

1.6fa73E+00 

l.e232E+CC 

3.6aC7E-02 

5.591AE-01 

1.5771E-02 

2.88266+00 

13 

1.O9C3E+00 

1.60756+00 

A. Vl37t-02 

9.57096-01 

1.71706-02 

2.83126+00 

lA 

l.b93t£+00 

1.5b72E+C0 

7,A7a56-J2 

9.5501E-01 

1.59566-02 

2.8727E+00 

15 

I.b9£6t+0C 

1.5896E+CC 

643C32E-02 

9.5321E-01 

1.2273E-02 

3.0295E+C0 

16 

1.6=P£iE+0C 

1.5327E+0C 

9*761AE-02 

9.52COE-01 

6. 66426-03 

3.0854E+00 

17 

1.69O2E+0C 

l.iClEE+OO 

l4l635b-01 

9.5157E-01 

l.Z517E-iO 

3.0244E+00 


Figure 26 — Continued 


95 




2-STATIGN no* ^0 Z« !?• 

4« 1 

3 69 




K 

Phi 

PPESSOP: 

lENSITT 

(.'-VEL 

V-VEL 
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APPENDIX 


CARD INPUT DATA FOR PROGRAM 4 

The input data for program 4 (called CHAOS, ref 4), are described 
below. The data produced by programs 1 and 2 and processed by program 5 
are read by a subroutine called FVDAT, while the data produced by pro- 
gram 3 are read by a subroutine called SCDAT Note that ENGLISH UNITS 
are used in these data. 


CARD INPUT DATA DEFINITIONS 


CARD INPUT DECK 


Card 

Format 

1 

(13) 

2 

(313) 

3 

(213) 

4 

(213) 

5 

(213) 

6 

(213) 

7 

(313) 

8 

(2F10. 3) 

9 

(2F10. 3) 

10 

(3F10. 3) 

11 

(3F10. 3) 

12^ 

(3F10.3) 

13 

(313) 

14 

(213) 

15 

(213) 

13A^ 

(313) 

14A 

(213) 

15A 

(213) 

i 

(A30) 

17 

(A30) 

18 

(13) 

19 

(8(13, 7X) ) 

20 

(13) 

21 

(8F10.4) 


Input Parameter 

JTOTAL 

NMTAPE, NSTAPE , DIMENSION 
NI, lOUT 
NPHIjNREC 
IN (5) , IN(6) 

JJR, NPSASF 

NPRINT , NPLOT , NPUNCH 

PSCALE, PHIMAX 

ALPHA, SCALE 

FS (13) ,FS (3) ,FS (2) 

FS (5) ,FS (6) ,FS (8) 

(PHIAO(J) ,CDEL(J) ,CEX(J) , J=1,NPHI) 
NMTAPE , NSTAPE , DIMENS ION 
NI, lOUT 
NPHIjNREC 

NMTAPE , NS TAPE , DIMENS ION 

NI, lOUT 

NPHI,NREC 


HEADl 

HEAD2 

NUM 

(J1(J) , J=1,NUM) 
NO 

(ZL(J) , J=1,N0) 


^PHI cards are read. 

^Cycle bach to card 13 as many times as required for 2 nrec = JTOTAL - 4 
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DATA 

Typical 

DEFINITIONS 

Name 

Value 

Description 

ALPHA 

30. 

FS(1) angle of attack (deg.) 

CDEL 

o 

o 

Cl = boundary layer thickness 
parameter for delta = C1*S**C2 

CEX 

0.0 

C2 = boundary layer thickness 
parameter for delta = Cl*S**C2 

DIMENSION 

17 

index giving array size for input 
data read in MCDAT ( = NUMK from 
CTI) 

FS (2) 

1716. 

gas constant 

PS (3) 

1.4 

gamma, specific heat ratio 

FS{5) 


PINF, free stream pressure, psf 

FS(6) 


RHOINF, free stream density, 
slug/ft^ 

FS (8) 


VINE, free stream velocity, fps 

FS (13) 

1.4 

gamma- INF, free stream specific 
heat ratio 

HEADl 

'SHUTTLE 147' 

plot heading 

HEAD 2 

'PERFECT GAS' 

plot heading 

IN{5) 

0 

gas type index, 0 for perfect, 
and -1 for equilibrium gas, and 

1 for nonequilibriiim (note that 
this IS different from programs 1 

2 and 3 ) 

IN (6) 

0 

gas file number for equilibrium 
gas, 2 for air 

lOUT 

2 

output control index for boundary 
layer data 

J1(J) 


array of station numbers for heat 


transfer crossplots, Jmax “ 
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DATA DEFINITIONS (Continued) 


Name 


Typical 

Value 


Description 


JJR 

JTOTAL 

NI 

NMTAPE 

NO 

NPHI 

NPLOT 


NPSASF 


NPUNCH 


NKEC 


0 number of points above the body 

surface used for edge conditions 

total number of points to be calcu- 
lated, four plus the number read 
from various tape files 

index to control type of input data, 
1 for SCDAT, 2 for MCDAT data, 
and 3 for PVDAT data 

nimber assigned ho the tape drive 
for input data from MCDAT 

total number of interpolated 
stations at constant span loca- 
tions, ZL 

19 number of streamlines 

index to control type of data 
plotted, 0 = no plot, 1 = print 
plot only, 2 = print plot, calcomp, 
3 = print plot, calcomp with dash 
lines, 4 = print plot, calcomp, 
data printout, 5 = contour plot 

0 to fill Newtonian streamlines to 
stag point, and 1 to start 
directly from input data 
\ 

0 for no cards punched and 1 
for punching cards for boundary 
layer input data 

number of records of data to be 
read in for the present file on 
tape drive NMTAPE 
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DATA DEFINITIONS (Concluded) 


Name 

NSTAPE 

NIM 

PHIAO 

PHIMAX 

NPRINT 

PSCALE 

SCALE 

ZL 


Typical 

Value Description 

number assigned to the tape drive 
for input from SCDAT or FVDAT 

number of values for Jl 

initial circumferential position 
of streamlines, deg. from leeward 

90. value of $ for terminating plots 

of streamlines 

0 for debug print, 1 for standard, 
and 2 for special print in CHAOS 

used for scaling of length dimen- 
sions for plotting 

used for scaling boundary layer 
data, not used now 

location for interpolated data at 
constant span, mnnber of values = NO 
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